Femtosecond pulses and dynamics of molecular photoexcitation: RbCs example 
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We investigate the dynamics of molecular photoexcitation by unchirped femtosecond laser pulses 
using RbCs as a model system. This study is motivated by a goal of optimizing a two-color scheme 
of transferring vibrationally-excited ultracold molecules to their absolute ground state. In this 
scheme the molecules are initially produced by photoassociation or magnetoassociation in bound 
vibrational levels close to the first dissociation threshold. We analyze here the first step of the 
two-color path as a function of pulse intensity from the low-field to the high-field regime. We use 
two different approaches, a global one, the 'Wavepacket' method, and a restricted one, the 'Level by 
Level' method where the number of vibrational levels is limited to a small subset. The comparison 
between the results of the two approaches allows one to gain qualitative insights into the complex 
dynamics of the high-field regime. In particular, we emphasize the non-trivial and important role 
of far-from-resonance levels which are adiabatically excited through 'vertical' transitions with a 
large Franck-Condon factor. We also point out spectacular excitation blockade due to the presence 
of a quasi-degenerate level in the lower electronic state. We conclude that selective transfer with 
femtosecond pulses is possible in the low-field regime only. Finally, we extend our single-pulse 
analysis and examine population transfer induced by coherent trains of low-intensity femtosecond 
pulses. 

PACS numbers: 33.80.-b, 34.80.Gs, 31.10.+Z, 33.15.-e, 



I. INTRODUCTION 

Rb and Cs atoms have been simultaneously trapped and 
laser cooled in a magneto-optic trap down to ultracold 
temperature (~ lOO^K). Ultracold RbCs molecules have 
been formed through photoassociation in excited vibra- 
tional levels of the Rb(5s)Cs(6p!/ 2 ) + , 0~ or 1 symme- 
tries. These molecules decay through spontaneous emis- 
sion, mainly toward stable levels of the Rb(5s)Cs(6s) 
a 3 S + electronic state; the upper of those levels has a 
binding energy in the range of 5 cm -1 [1]. The relevant 
molecular terms are shown in Fig. 1. 
In the heteronuclear RbCs molecule, two-step conver- 
sion processes from the a 3 S + state (denoted below by 
'a') toward the X 1 E + state (denoted below by 'X') are 
possible by using, as intermediate step, levels of the + 
or 1 symmetries, with a spin-mixed character. As a re- 
sult, molecules in the absolute ground level Rb(5s)Cs(6s) 
X 1 Y 1 + v" = are formed. These processes have been re- 
cently investigated experimentally [1, 2] and theoretically 
[3-5]. 

Ultracold stable polar molecules in their absolute ground 
vibrational level have been populated for the first time 
[1, 2] using a two-color incoherent population trans- 
fer through a low-lying level of the 1 state. A reso- 
nant 'pump' laser pulse transfers the population of the 
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metastable, vibrationally excited a 3 S + molecules to an 
electronically excited level; then a second tunable 'dump' 
laser pulse resonantly drives the population to the ab- 
solute ground level. The two laser pulses used in this 
stimulated transfer have a duration of about 5 ns. 
In the KRb molecule, using a Stimulated Raman Adi- 
abatic Passage (STIRAP) with counterintuitive pulses 
in the microsecond range, Ni et al [6] transferred ex- 
tremely weakly-bound Feshbach molecules in the a elec- 
tronic state toward the lowest vibrational level either of 
the stable X or of the metastable a states using interme- 
diate level with symmetry 1. 

For several years, researchers at Aime Cotton Labora- 
tory are exploring theoretically, on the example of the 
Cs2 and Rb2 molecules, coherent schemes using chirped 
laser pulses to form molecules in an excited electronic 
state through photoassociation of ultracold atoms, and 
then to stabilize them through stimulated emission [7-9] . 
The motivation was to fully exploit optical techniques for 
controlling the formation of cold molecules in the abso- 
lute ground level. The studied laser pulses were in the 
picosecond range, the domain well-adapted to the vibra- 
tional dynamics of the wavepackets created by the pulse 
in the light-coupled electronic states. 
However, from a technological point of view, picosecond 
lasers and corresponding pulse shapers are not yet avail- 
able. On the other hand, in the femtosecond domain 
there were important recent developments of efficient 
laser sources and pulse shapers. Furthermore, coherent 
trains of pulses, obtained from mode-locked femtosecond 
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FIG. 1: (Color online) Photoexcitation of RbCs. Lower 
panel: Diabatic potentials for the stable .X" 1 E + (lower thick 
dashed black line), the metastable a 3 E + (lower thick contin- 
uous black line) and the excited fc 3 II (upper thick continuous 
black line) and ^4 1 E + (upper thick dashed black line) elec- 
tronic states. The origin of energy is set at the dissociation 
limit Rb(5s)Cs(6s). Vibrational wavefunctions for the ini- 
tial level a 3 E + v" = 37 (upper continuous yellow [clear-gray] 
line) and for the final level X l Y, + v" — (lower continuous 
yellow [clear-gray] line) in the two-color process. Wavefunc- 
tions for the fe 3 II v' = 43 level (continuous orange [dark-gray] 
line) corresponding to the resonant transition ('oblique' black 
continuous line) and for the fo 3 il v' — 195 level (continuous 
red [black] line) corresponding to the off-resonant 'vertical' 
transition (vertical dashed black line) . The wavefunctions are 
drawn at their absolute energy. The bandwidth of the laser 
pulse is shown as the thick black horizontal line. Upper panel: 
Adiabatic potentials + P1/2 (thick continuous black line) and 
+ P3/2 (thick dashed black line) . The energy origin is at the 
Rb(5s)Cs(6p) dissociation limit. The two vibrational com- 
ponents of the coupled + v' = 70 level in the Hund's case 
c representation are drawn, respectively fe 3 n (medium-thick 
red [black] line) and J 4 1 E + (thin green [medium-gray] line). In 
the inset, these two components are compared with the wave- 
function of the fe 3 II v' = 43 level (medium orange [dark-gray] 
line) for 8 < R < 12 a . 



lasers [10], permit a transient coherent accumulation of 
population, manifested by the enhancement of transition 
probabilities and by a gain in the spectral resolution [11]. 
Our objective here is to analyze the possibilities offered 
by femtosecond sources in implementing efficient two- 
color paths for transferring vibrationally-excited ultra- 
cold molecules to their absolute ground state. In this 



scheme the molecules are initially produced by photoas- 
sociation or magnetoassociation in bound vibrational lev- 
els close to the first dissociation threshold. Numerical 
analysis is carried out for the RbCs molecule. More pre- 
cisely, the present paper is devoted to the choice of the 
optimal pulse for implementing the first step of the two- 
color paths. Notice that femtosecond pulses have a broad 
bandwidth and may reach high intensities. Consequently 
we have to analyze the dynamics of coherent excitation 
of a large number of vibrational levels, from the low field 
up to the high field regime. 

To solve the time-dependent Schrodinger equation, we 
first use the 'Wavepacket' method (WP), where we cal- 
culate globally the evolution of vibrational wavepackets 
propagating along electronic states coupled by the laser 
pulse [7]. Using this approach, it appears that, in the 
high-field regime, the calculated dynamics and the popu- 
lation transfer drastically differ from what is expected 
from intuitive two-level-system arguments. To under- 
stand these surprising results, we compare the WP re- 
sults to solutions obtained using a small subset of vi- 
brational levels: we refer to this model as the 'Level by 
Level' method (LbyL). In both approaches, the depen- 
dence of the wave function on the interatomic distance R 
is obtained from the Mapped Fourier Grid Hamiltonian 
(MFGH) method [12, 13]. 

By comparing the WP results with the LbyL solutions, 
we precisely identify vibrational levels critically responsi- 
ble for the strongly nonlinear dynamics in the high-field 
regime. In the high- field regime, the dynamics of the pho- 
toexcitation process is governed both by nearly- resonant 
and by far-from-resonance excitations. The adiabaticy 
of the resonant and non-resonant excitations can be eas- 
ily analyzed in detail in the simple case of a two-level 
system. For a multilevel system, we show that, in the 
high-field regime, the dynamics of time-evolution of the 
population in nearly-resonant levels is strongly affected 
by the adiabatic excitation of far-from-resonance levels. 
For a particular level, the adiabaticity of the excitation 
by an unshaped Gaussian pulse is found to be simply 
related to the value of its detuning with respect to the 
carrier laser frequency. In the photoexcitation process 
under study, the initial level lies close to the dissocia- 
tion threshold, in an energy domain where the density 
of vibrational levels is high. We show that the presence 
of such a quasi-degenerate group of levels in the ground 
electronic state leads in the high field regime to a spec- 
tacular blockade of the excitation process. 
We conclude from the analysis that, while femtosecond 
laser pulses are concerned, control of the photoexcitation 
process is possible only in the low field regime. To im- 
prove the efficiency of the population transfer, we investi- 
gate some schemes using coherent trains of low-intensity 
femtosecond pulses. 

The paper is organized as follows. First we specify the 
photoexcitation process (Sec. II A) and also characterize 
Gaussian pulse (Sec. II B). Then we briefly describe the 
two employed approaches (the WP and LbyL methods) 



to solving the time-dependent Schrodinger equation (Sec. 
IIC). The photoexcitation dynamics is dramatically af- 
fected as the pulse intensity is increased. It's dependence 
on the pulse intensity is computed in the WP approach 
and is described in Sec. III. These results are further an- 
alyzed in Sec. IV in the framework of the LbyL method. 
This framework allows us to identify levels responsible 
for the observed photoexcitation dynamics (Sec. IV A). 
We further exhibit the link between adiabaticity and de- 
tuning first in the simple case of a two-level system (Sec. 
IV B) and then for the multi-level system under study 
(Sec. IV C). The excitation blockade due to the presence 
of quasi-degenerate group of levels in the ground state 
is studied in Section IV D. Finally, we comment on the 
photoexcitation dynamics induced by coherent trains of 
low-intensity femtosecond pulses in Section V. 

The paper contains several appendices used for recapit- 
ulating essential results and to precise notation. Ap- 
pendix A briefly reviews the Mapped Fourier Grid Hamil- 
tonian (MFGH) employed throughout the paper. The 
'Wavepacket' and the 'Level by Level' methods are de- 
scribed in the Appendix B. Appendix C recalls the def- 
inition of the diabatic and adiabatic bases used in our 
analysis. A simple model for the blockade of excitation 
due to the presence of a quasi-degenerate group of lev- 
els in the lower electronic state is described in Appendix 
D, whereas Appendix E lists relevant properties of ultra- 
short pulse trains. 



II. PHOTOEXCITATION OF RbCs 
A. Photoexcitation process 

In the RbCs molecule, it has been shown that the two- 
color path a 3 £+w" = 37 -> 0+v' = 70 -> X^+v" = is 
very efficient in transferring to the absolute ground level 
AT 1 S + f " = the molecules obtained in the av" = 37 level 
after photoassociation followed by spontaneous radiative 
decay [14]. The + symmetry results from the coupling 
through the spin-orbit interaction of the singlet A 1 'E + 
and the triplet b 3 H electronic states. The + v' = 70 
level is a mix of vibrational levels b 3 ILv' (52.7%) with 
v' ~ 43 and of A^+v' levels (47.3%) with v' ~ 25. 
In the first step of the two-color path, only the |6 3 IIu') 
components of the coupled wave functions |0 + u') can be 
excited; we have shown that the excitation probabilities 
av" = 37 ->• W = 43 and av" = 37 -> 0+v' = 70 level 
are very similar. Therefore, in this paper, we restrict the 
analysis of the photoexcitation dynamics to the study of 
the a 3 £ + u" — )• b 3 Hv' transition. The rotational structure 
of the vibrational levels as well as the hyperfine structure 
are ignored. 

We consider excitation by a Gaussian laser pulse with a 
duration tl and a carrier frequency ujl/2it resonant with 
the transition between the vibrational levels av"o = 37 



and b 3 Uv' = 43, 

hu L = E(b 3 nv' a = 43) - £(a 3 £+u" = 37), 



(1) 



where £(a 3 £+v" = 37) and E(b 3 Iiv' a = 43) are absolute 
energies of the two levels. 

The initial level has a binding energy of only 5.52 cm -1 
and it lies very close to the Rb(5s)Cs(6s) dissociation 
limit. The excited level with binding energy 4392 cm -1 
with respect to the Rb(5s)Cs(6p) dissociation limit is 
tightly bound (Fig. 1). There are substantial differences 
in the two vibrational wave functions. The wave func- 
tion of the initial level a 3 £ + vo" == 37 extends from 9 
to 27 ao (ao denotes the Bohr radius) and the wave 
function of the resonant level b 3 IL v' Q — 43 is located 
at much smaller internuclear distance, 7 ao to 11 ao. 
As a result, the Franck-Condon factor is relatively small 
(| (au " = 37|6v = 43) | 2 = 1.16 x 10~ 3 ). 
In the same Fig. 1 we also show the wave function, in 
the Hund's case a representation, of the spin-orbit-mixed 
vibrational level + v' — 70, which has an energy close 
to the energy of the pure Hund's case a resonant level 
bv' = 43. One should notice the similarity between the 
vibrational component in the b triplet state of the wave 
function + v' = 70 and the vibrational wave function of 
the pure bv' — 43 level for 9.5ao < R < llao, that is in 
the i?-range where the overlap of both wave functions is 
the largest. 

The wave function of the bv' = 195 level, strongly off- 
resonant with the studied laser pulse but connected to 
the a «o" level through a 'vertical' transition (the outer 
turning points of both wave functions are located at 
Rout ~ 26 ao), is also reported in Fig. 1. The correspond- 
ing Franck-Condon overlap, | (avo" = 37\bv' — 195) | 2 = 
0.183, is much larger than that one of the resonant tran- 
sition. 



B. Characteristics of the laser pulse 

The laser pulse is assumed to have a Gaussian profile and 
to be Fourier-transform-limitcd, with a time-independent 
carrier frequency fixed to wj,. We do not consider chirped 
pulses because the mechanism of adiabatic population 
transfer occurring during excitation with chirped pulses 
has been previously extensively analyzed and optimized 
[7, 8, 15, 16]. The motivation of the present work is to 
investigate a completely different excitation mechanism, 
resulting from the use of ultrashort unchirped pulses, and 
to interpret in detail its dynamics. 

The laser pulse is described by an electric field with an 
amplitude £{t) varying with time as: 

£(t) = £ f(t)cos[u L t] = £(t)+£*(t) 

- ^/(t)exp[*W£t] + y/WexpHwit] , (2) 

where £q is the maximum amplitude and £_(t) denotes 



the complex time-dependent amplitude. The Gaussian 
envelope /(£), with maximum f(tp) = 1, is given by 



f(t) = cxp 
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The instantaneous intensity I(t) of this pulse illuminat- 
ing an area a, is equal to 



TU , ^puisc /4m 2 
I{t) = -z^\l^^exp 
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where I(tp) = II = ceq£q/2 (c is the velocity of light, eo 
the vacuum permittivity). I(t) has a full width at half 
maximum (FWHM) equal to tl ■ The pulse duration and 
the energy -E pu isc of the pulse satisfy: 
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In the spectral domain, the electric field £(u>) is ob- 
tained from the Fourier transform of the complex time- 
dependent electric field £_(t), 
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(6) 



For the pulse of duration <~ 100 fs considered here, the 
bandwith Su = 41n2/r L , defined by the FWHM of 

\£(u> — ojl)\ , is of the order of ~ 150 cm -1 . 



C. Photoexcitation dynamics: 'Wavepacket' and 
'Level by Level' descriptions 

To analyze the dynamics of the photoexcitation process 
(Eq. (1)), we consider the time-dependent Schrodingcr 
equation describing the internuclear dynamics of the Rb 
and Cs atoms 

[H rnol -j2-E(t)]*(t)=ih^-*(t), (7) 

at 

where H mo i denotes the molecular Hamiltonian in the 
Born-Oppenheimer approximation and where the cou- 
pling between the laser and the molecule, written in the 
dipole approximation, is expressed in terms of the dipolc 
moment operator jl. The electric field of the laser pulse 
with polarization e p reads E(t) = £(t)e p . 
In the excitation process, we focus on the redistribution 
of the population between the vibrational levels, disre- 
garding rotational components of the wavepackcts ^(t). 



This approximation is justified because the centrifugal 
energy is negligible and thereby vibrational wavepackcts 
do not depend on value of the total angular momentum 
J. All our calculations were carried out for a fixed value 
of J, J = and below we do not identify it explicitly. 

In a simple model restricted to the ground g and excited 
e electronic states, the two radial components x 9 (i?, i) 
and Xe(-R)i) of the wavepacket ^(R, t) arc solutions of 
the coupled system: 



m dt\Xe(R,t) 



(8) 






V g {R) W ge {t)cos[uj L t] 
V e (R) 
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Xg(R,t) 
Xe(R,t) 



where V g (R) and V e (R) denote the potentials in the 
ground and excited states. The coupling of the two elec- 
tronic states can be written in terms of W ge (R, t) 

W ge (R,t) = -D ge {R)£ f{t) = -2W L f(t), (9) 

where D ge (R) denotes the electronic dipole transition 
moment resulting from the integration of ft over the elec- 
tronic wave functions of the ground and excited electronic 
states. We disregard the i?-dependence of the electronic 
transition dipole, which is taken equal to its asymptotic 
value D. Finally, Wl = D£o/2, the maximum strength 
of the coupling, is proportional to the square root of the 
maximum intensity 1^. 

The radial part of the wavepackets Xg(Rjt) (resp. 
Xe(R, t)) is a coherent superposition of the stationary 
vibrational wave functions, eigenstates <p giV » (R) with 
energy E g y (resp. ip e ^ v /(R) and E e , v i) of the time- 
independent Schrodingcr equation involving the poten- 
tial V g (R) (resp. V e (R)). Numerically the radial depen- 
dences of all functions are described by using the Mapped 
Fourier Grid Method (MFGH) [12, 13]. Let us empha- 
sized that for a single potential, the eigenstates consist 
of bound levels and discretized scattering levels, which 
are automatically included in the decomposition of the 
wavepacket (see Appendix A). A spatial grid of length L 
with TV mesh points is used for each potential yielding a 
quasi-complete set of N cigenfunctions (see Ref. [17]). 

Two methods are used to solve the time-dependent 
Schrodingcr equation in the rotating wave approxima- 
tion (RWA). The first method, the Wavepacket descrip- 
tion, consists in determining directly the vibrational 
wavepackets Xg {R> t) an d Xe (R, t) created by the laser 
pulse on both electronic states g and e. Studying the ex- 
citation from the vibrational level a 3 X + «o" = 37 (Fig. 1), 
the initial state is chosen to be this initial vibrational 
level: X g (.R, t = 0) = f a ,v ( r and Xe{R,t = 0) = 0. De- 
tails on the numerical methods, presented in Refs [8, 9], 
are summarized in Appendix Bl. The time-dependent 
Schrodingcr equation is solved by expanding the evolu- 
tion operator in Chebyschev polynomials [18]. With the 
MFGH method being used to represent the radial de- 



pendence of the wavepackets, the WP method is a global 
approach which automatically incorporates contributions 
of the complete set of vibrational levels (p g , v >< (R) and 
ipe.y (R) with < v', v" < N - 1 . 

The second approach, the Level by Level description, an- 
alyzes the coupling by the laser pulse of some beforehand 
selected subsets of vibrational levels g , e m , with g n and 
e m being numbers of levels in the ground and excited 
state vibrational subsets respectively. 

The (g n + e m ) chosen levels result in the formation of the 
ground and excited wavepackets written, in the 'interac- 
tion representation' [19], as: 
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where the phase factor accounts for the 'free evolution' 
of the stationary vibrational levels. In the RWA approx- 
imation, the instantaneous probability amplitudes a v >< (t) 
and b v , (t) are determined by solving a system of {g n +e m ) 
coupled first-order differential equations (Eq. (B6)) pre- 
sented in the Appendix B2. For the initial state of the 
system, the probability amplitude of the a 3 £ + «o" = 37 
level is set to unity: a v » (t = 0) = 5(v" , w " ) and 
b v ,(t = 0) = for all the considered v' values. The rel- 
evant molecular structure data are the relative energies 
A^„ Vo „ for the ground levels (resp. A£, v , for the excited 
levels) with respect to the resonant level vq" (resp. v' ), 



A 9 ,, „ 
A e , , 



— &g,v 

= E, 
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and the overlap integrals 



(v'\v") = (<p e ,v>\<Pg,v")- 



(11) 



(12) 



Notice that because of the resonance condition, the en- 
ergy spacings A^„ v „ and A®, v , may be expressed in 
terms of the detunings 5 v ' tV » (Eq. (B7)). 

The WP and LbyL methods are compared in Appendix 
B 3. The WP/MFGH approach allows one to expand the 
wavepackets \g {R> t) and Xe (R, t) over the complete set 
of N vibrational levels of the g and e electronic states: 



Xg(R,t) 



Xe(R,t) 



E„ 



y^ a v » (t) exp[-i 9 ' v t] ipg. v " (R) 



=0,N-1 



E„ 



y^ b v i(t)exp[-i—Y-t]<p e , v i(R).(lS) 



=0,N-1 



P e (t) = ( Xe (R,t)\Xe(R,t)) , 

Pg(t) = (Xg(R,t)\ X g(R,t)) . 



(14) 



More detailed information is provided by decomposing 
the wavepackets in the basis of unperturbed vibrational 
levels v' or v" of both electronic states e org, 

Pgv(t) = K^,»-(^)IX fl (^*)>| 2 = |Ot;-(*)| 2 . 
Pev'(t) = \(p e y(R)\Xe(R,t))\ 2 = K(t)\ 2 , (15) 

which gives the instantaneous population of each station- 
ary vibrational level. 

For the LbyL approach, populations similar to those 
defined in Eqs. (14) and (15) can be introduced. 
Naturally the LbyL approach is equivalent to the WP 
description if and only if the sets g and e m encompass 
complete sets with g n = e m = N levels, that is all bound 
levels and all levels of the discretized continua (Appendix 
A). We emphasize that the WP description automati- 
cally takes advantage of the completeness of the set of 
cigenfunctions provided by the spatial representation of 
the Hamiltonian on a grid. Furthermore, the description 
of the dynamics does not depend on the choice of the 
grid, provided that a sufficiently wide domain of energy 
is covered by the eigenvalues obtained in the MFGH di- 
agonalization of the Hamiltonian matrix. 



III. WAVE PACKET DESCRIPTION: FROM 
LOW FIELD TOWARD tt-PULSE 

A. 7r-pulse condition 

Our goal is to find a pulse which yields a population 
transfer as large as possible from the initially popu- 
lated vibrational level a 3 £ + vq" — 37 toward the vi- 
brational level 6 3 n v' = 43. As mentioned above, 
we consider only the case of an unchirped transform- 
limited Gaussian pulse, resonant with the transition 
a 3 £ + i>o" = 37 — >• 6 3 1Ti>q = 43, with a duration in the 
femtosecond domain. The chosen duration is r^ = 120 



The evolution of the total population in the two elec- 



fs, much smaller than the vibrational period T!"; 6 „ = 
Airh/(Eg^ Vo -> + i) — -E g ,(t, »-i)) ~ 22 ps for the initial level 
a 3 £ + i>o" = 37. It is only 6 times smaller than the vi- 
brational period T mb , = 0.72 ps in the excited state. 
Consequently, in the excited electronic state, there are 
only 6 nearly-resonant levels lying within the bandwidth 
6u> = 122cm~ 1 = 5.59 x 10~ 4 au of the pulse, the levels 
41 < v' < 45 with detuning S v ' Vo " respectively equal to 
-92.0, -46.0, 0, +45.8, +91.5 cm" 1 . 

The pulse is characterized by the electric field amplitude 
£q or, cquivalently, by the pulse intensity II or by the pa- 
rameter Wl (Eq. (9)). Given a pair of levels (say i>o" and 
v' ) we may also introduce the accumulate pulse area [20] 



e(t) = D£ |KK')| / f(t')dt', 



(16) 



where (v' \vo") denotes the overlap integral of the res- 
onant transition (Eq. (12)). The total pulse area of a 
Gaussian pulse is 

Q F = 0(+oo) = W L T L \(v' \v »}\^/ —. 

In a two-level system, the angle 8(t) fully determines the 
probability amplitudes of the lower level a(t) and of the 
resonantly-excited (i.e., when S v > „ » ~ ) level b(t) as 

[20] 
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, 6(i) = i sin 
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(17) 



The 7r-pulse for a resonantly-driven two-level system is 
defined as 0_f = w, 



tlW l Hk")l = \M^ ov 4*/— fi(t P ) = Su,(l8) 



where £l(tp) = Wl \{vq\vo"}\ is the Rabi coupling (see 
Eq. (20)) for the resonant transition v' — > «o" at the 
pulse maximum t = tp . 

Accounting for the overlap integral |(uol u o"}| = 
0.03462 au. and for the pulse duration tl = 0.12 ps= 
4961.11 au, the 7r-pulse condition is satisfied when 



Wl = -D£ = 6.076 x 10~ 3 au or I L = 493GW/cm 2 . 

This large value of intensity is due to the small value of 
the overlap integral and to the short pulse duration. 
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FIG. 2: (Color online) WP approach: variation with time 
(in ps) of the population in the excited electronic state 6 3 II. 
Panels a), b), c): low field excitation Wl = W£/120. Panels 
d), e), f): high field excitation Wl = W£. Panels a) and 
d): total population ^„/Z |6„'(t)| 2 in the excited electronic 
state. Panels b) and e): population |&„/ =43 (t)| 2 in the reso- 
nant level Dq = 43. Panels c) and f): population |b„/ =1 g 5 (£)| 2 
in the far from resonance level v' — 195, corresponding to the 
'vertical transition' defined in Fig. 1. Duration of the pulse 
[tp — TL,tp + tl] (vertical dashed blue [black] lines) maxi- 
mum at t = tp (vertical thick continuous blue [black] line). 
The populations have been multiplied by the factor indicated 
in the upper right corner. 



the initial level is almost not modified during the pulse. 
After the end of the pulse, for t 3> tp+r L , the population 
of the level v' in the excited electronic state is equal to: 



B. Low field excitation 



\b v ,{t->+oo)\' 



^-(v'\vo"Y 



We first consider a weak pulse, Wl = 5. 10~ 5 au, with 
a pulse area @p — 7r/120, corresponding to an intensity 
at the maximum of the pulse II — I(tp) = 34MW/cm 2 . 
The initial population in the a 3 £ + i>o" = 37 level is set 
equal to unity. The evolution with time of the total pop- 
ulation in the excited electronic state & 3 II and in the res- 
onant level v' = 43 is reported in Fig. 2a, b. The consid- 
ered populations increase monotonously during the pulse 
and the total transfer is very small (0.000343), with half 
population (0.000161) in the resonant level v' = 43. For 
the v' = 42 and v' — 44 levels, which have a detun- 
ing with respect to the central laser frequency smaller 
than 5oj/2, the population at the end of the pulse is re- 
spectively 0.000076 and 0.000088. There is almost no 
population in the levels v' < 40 or v' > 46. 
In the perturbative limit, the amplitude of population of 
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where S v > yV0 » is the detuning of the excitation of the ev' 

level from the gvo" level and where \£(u> — ojl)\ 2 (Eq. (6)) 

is the spectral density of the pulse. 

In this limit, the population transferred from the level 

gvo" toward the level ev' is proportional to the Franck- 

Condon factor (v'\v Q v ) and to spectral density of the 

pulse at the excitation frequency [21]. As a result, for 

the weak perturbative pulses, only the nearly-resonant 

levels, such as |5 O ' )U0 "| < 5u>, are excited. 

The population distribution in the vibrational levels is 

presented in the left column of Fig. 3 for the excited elec- 
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FIG. 3: (Color online) WP approach: population |fej,/(i)| 2 in 
the levels 6 3 HV (panels a) and c)) and population \a v » (t)\ 2 in 
the levels a 3 E + v" (panels b) and d)), as a function of the en- 
ergy of the corresponding levels (in cm -1 ) with respect to the 
Rb(5s)Cs(6p) and Rb(5s)Cs(6s) dissociation limits respec- 
tively. The distributions of population are shown either at 
the maximum of the pulse t — tp = 0.6 ps (black dots) and 
after the end of the pulse t = 2 ps (orange [gray] squares). 
Panels a) b): low field excitation Wl = Wl/120. Panels c) 
d): high field excitation Wl — W£. The resonant level v'q is 
indicated by a large full black circle at t = 0.6 ps and a large 
full orange [gray] square at t = 2 ps, the level bv' — 195 is 
represented by a full black triangle. 
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FIG. 4: (Color online) WP approach. Panel a): variation of 
the population \b v /(t — > +oo)| 2 remaining after the pulse in 
the nearly-resonant excited levels & 3 HV, as function of the 
laser coupling Wl (in units of 10~ 3 au); v' = 41 (medium- 
thick red [black] line); v' = 42 (thin black line); v = v'q = 43 
(thick black line); v' — 44 (medium-thick clear-gray line); 
v' = 45 (medium-thick orange [gray] line). The low field 
excitation WX/100 is indicated by the thin dashed vertical 
red [black] line. The couplings Wl x — 4.5 x 10~ 4 au (black 
line), Wl 2 = 1-15 x 10 -3 au (orange [gray] line) and Wl 3 = 
1.75 x 10 -3 au (clear-gray line) corresponding to the maxima 
of |b„/ =43 (i — ► +oo)| 2 are indicated by the vertical arrows 
at the top of the panel. Right panels: couplings WL t i = 
1 — 3. Panel b): population transferred to the levels b s Tlv' 
as function of v' . Panel c): Population redistributed in the 



levels a 3 E 



3y+„ 



as function of v" . 



tronic state (panel a) and for the lowest electronic state 
(panel b), either at the maximum of the pulse (t = tp) 
or after the end of the pulse (t — 2 ps). The population 
of the excited vibrational levels v' ~ 180 — 200 always 
remains smaller than that of the nearly-resonant levels 
40 < v' < 45, and, at the end of the pulse, only these 
levels remain populated. In the low-field limit, the dy- 
namics of the excitation process involves almost only the 
nearly- resonant levels (Figs. (2) and (3)). 



C. Increasing the field strength 

Now we vary the laser coupling Wl, and explore the pop- 
ulation \b v /(t — > +oo)| 2 transferred to the excited levels 
bv' with 41 < v' < 45. The results of our WP cal- 
culations are shown in Fig. 4a. In the low-field limit, 
the populations increase proportionally to W^, and, as 
already noticed, only the levels v' = 42, 43 and 44 
are significantly populated. However, when the pulse 
area/intensity are increased, the population in the lev- 
els with v' < 41 or v' > 45 becomes comparable to 
the population in the nearly-resonant levels. The pop- 
ulation in the resonant level at the end of the pulse, 



\b v > (t — > +oo)| 2 , first increases with increasing Wl and 
reaches, for W L - WJ/14 = 0.000425 au, a relatively 
small maximum, 0.0052 ~ 1/14 2 . This coupling corre- 
sponds for the resonant transition to an 'effective' pulse 
area of 7r/14, still in the low-field regime. As Wl is in- 
creased further, the \b v > (t = +oo)| 2 oscillates with a pe- 
riod roughly equal to AWl = 0.0007 au. Notice that as 
a function of Wl, the values of the population maxima 
decrease after two oscillations. This behavior strongly 
differs from what one would expect intuitively for the 
resonantly-excited two-level system [<?t>o" , e ^o] : i n that 
case, the population would oscillate between the values 
of and 1, with a period equal to 2W[ : the value of 1 



Wl 



6.076 x 10"^ au. 



being reached at Wl 

The population distribution among the levels of the ex- 
cited 6 3 II and initial a 3 £ + electronic states after the pulse 
is presented in Fig. 4b, c for three values of the coupling 
Wl- These couplings correspond to the first three max- 
ima in the variation of |6„' = 43(£ — > +oo)| 2 as a function 
of Wl (see the vertical arrows at the top of Fig. 4a) . For 
Wl x — 4.5 x 10~ 4 au, only three nearly-resonant levels 
are populated and no significant redistribution of popu- 



lation occurs in the av" levels. For Wl 2 = 1.15 x 10~ 3 
au, more bv' levels, with 39 < v' < 46, are populated and 
the population is recycled back to levels av" of the initial 
state with 25 < v" < 43. For W Ls = 1.75 x 10~ 3 au, a 
still larger number of au" and bv' levels is involved in the 
redistribution of population. 



D. 7r pulse: resonant and far-from-resonance 
excitation 

The time-evolution of the total population ^2 V ,Z \b v > (t)\ 2 
transferred to the excited electronic state 6 3 II during the 
excitation by a pulse with a large coupling strength W[ 
is presented in Fig. 2d. Population maximum (0.094) 
is attained at the maximum of the pulse t = ip; it be- 
comes smaller when the pulse intensity decreases. The 
final value, equal to 0.019, is much smaller than unity. 
The evolution of the population \b v i (£)| 2 of the resonant 
level v' = 43 is shown in Fig. 2e. This population docs 
not increase monotonically, as one would expect for a tt- 
pulse in a two-level system, but exhibits several (~ 11.5) 
oscillations and the transfer is low (0.00064). A sim- 
ilar behavior is observed for the nearly-resonant levels 
v' = 42 and v' = 44 with final populations of 0.00043 
and 0.00059, respectively. Figure 3 shows the popula- 
tion distribution over various levels of the excited (panel 
c) and of the lowest (panel d) electronic states at two 
times t = tp = 0.6 ps and at t = 2 ps. We find that at 
the end of the pulse a significant fraction of the popula- 
tion is transferred to a large number of strongly-bound 
b 3 H v' levels, mainly to the levels 26 < v' < 56 with bind- 
ing energies in the range of -5200 to -3800 cm -1 . The 
most populated levels, v' ~ 31 and v' ~ 51, with respec- 
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cm -1 , have a population ~ 0.0013, equal to twice the 
population of the resonant level v' = 43. Population 
is also redistributed within bound and scattering levels 
of the ground a 3 £ + v" state, in particular within levels 
32 < v" < 41 (population > 0.005). The difference in 
the energies of these levels with respect to the initially 
populated level A^„ Va „ = 8 v i iV « (Eq. (11)) lies in the 
range —4.2 cm -1 < S v > „» < 11.3 cm -1 . Only 76% of the 
population remains in the initial w " = 37 level. 
At the maximum of the pulse, there are many levels of 
the excited electronic state, v' ~ 180 — 200, which have 
a population larger by a factor of at least 10 than the 
population in the nearly-resonant 41 < v' < 45 levels 
(\8 V ',v " | < 100 cm -1 ). These strongly-populated levels 
are such as 5 v i. Vo » ~ 4300 cm -1 , so they lie far outside 
the pulse bandwidth and correspond to highly-far-from- 
rcsonance excitations. Because of their high population 
during the pulse, these levels contribute significantly to 
the excitation dynamics. The time evolution of the pop- 
ulation |6</ = i95(£)| 2 of the b 3 Hv' — 195 level, is reported 
in Fig. 2f. This is the most populated level in the excited 
electronic potential with a population reaching 0.0134 at 
the maximum of the pulse. The time-dependence of this 



population follows that of the envelope of the pulse in- 
tensity, |/(t)| 2 (Eq. (3)). Note that in the low-field case 
(Fig. 2c), the population of the v' = 195 level is always 
negligible (< 10~ 6 ). 

It is to be emphasized that this behavior can not be ex- 
plained as Rabi cycling, contrarily to what could be in- 
tuitively expected considering the large value of the in- 
stantaneous Rabi coupling £t(v', Vo" ,t) arising from the 
large value of the detuning 8 v ' = ig5 yVo » . We recall the def- 
inition of the instantaneous Rabi coupling at time t for 
non- resonant transition av" — > bv': 

Q(v', v" , t) = ^\S v , y , P + 4|/(i) (v>\v") W L \\ (20) 

where (v'\v") denotes the Franck-Condon factor. The 
b 3 Hv' = 195 level is the level of the excited electronic 
state possessing the largest population at the maximum 
of the pulse. This can be understood by reminding that 
this level is excited from the initial level a^o" through a 
vertical transition (see Sec. II A). 

The importance, in the strong field regime, of off- 
resonant excitation of levels strongly favored by high 
Franck-Condon factors but lying energetically above the 
spectral bandwith of the pulse has been experimentally 
observed in the photoassociation of ultracold atoms with 
shaped femtosecond pulses [22, 23]. 



IV. ANALYSIS OF THE tt-PULSE DYNAMICS: 
LEVEL BY LEVEL DESCRIPTION 

The WP results demonstrate that, for the high coupling 
strength W£ , the dynamics of the excitation process in- 
volves a large number of vibrational levels, both in the 
ground a 3 S + and in the excited 6 3 II electronic states. 
To better understand the dynamics of the population of 
these levels, we performed LbyL calculations, with var- 
ious subsets g and e m of bound and quasi-continuum 
(scattering) levels. These subsets are simply denoted as: 
[v 2 ",...v g :\v' 1 ,v' 2 ,...v' e J. 



A. Levels involved in the dynamics 

1. LbyL basis set reproducing the WP dynamics 

In the first step we try to reproduce, by optimizing the 
restricted LbyL basis set, the time-evolution of the total 
population transferred to the excited electronic state by 
7r-pulse (Wl = W£). Some representative results are 
displayed in the left column of Fig. 5, where the following 
basis sets are considered: set A: [v" — 30 — 50, v' = 
- 218], set B: [v v = 20 - 50, v' = 30 - 50, 190 - 200] 
and set C: [v" = - 204, v' = - 218]. These levels are 
either bound or discretized scattering vibrational levels 
in the a 3 £ + or 6 3 II electronic states. Let us remark, that, 
with the mesh grid used in the MFGH approach, only a 
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FIG. 5: (Color online) LbyL description with three different 
basis sets of levels in the lower (y") or the excited (v') elec- 
tronic states; these sets, labeled A, B and C, are defined in 
the text: variation of the population in the excited b 3 HV 
levels as a function of time (in ps). The pulse (vertical con- 
tinuous and dashed blue [black] lines) satisfies the 7r-pulsc 
condition. Panel a): total population for set A (continuous 
orange [clear-gray] line) , B (dashed red [dark-gray] line) and 
C (continuous black line). Panels b), c) and d): 'optimal' 
basis set C [v" = - 204, v' — — 218] reproducing the re- 
sults of the WP method (Fig. 2); b): total population in the 
bound levels v' = — 218; c): population in the resonant 
level v' = 43; d): population in the off-resonant level 6 3 IT 
v = 195, corresponding to the vertical transition. Some pop- 
ulations have been multiplied by the factor indicated in the 
upper right corner. 



small energy range (0 < E < 0.01 cm -1 ) is described by 
'physical' scattering levels (see Appendix A). 
The relatively large set B includes, in the lower state, 
bound levels lying close to the initial one, «o" = 37, and, 
in the excited state, levels located in the vicinity of the 
resonantly excited v' = 43 level or in the vicinity of the 
far-from-resonance v' = 195 level corresponding to the 
vertical transition. For this set, the total population at 
the maximum of the pulse t = tp is larger by a factor 2 
than the population obtained by using the WP approach. 
At the end of the pulse, a too large population (~ 0.12) 
remains in the excited state. For the set A, which in- 
cludes all the bound levels in the excited state and, in 
the lower state, a smaller number of levels located in the 
vicinity of the initially populated one, a similar behavior 
is obtained, yielding the same final population transfer, 
but a slightly smaller maximum value at t ~ tp. 
To reproduce in the LbyL approach the results obtained 
in the WP approach, we have found that it is neces- 
sary to employ the set C which includes all bound vibra- 



tional levels in the excited state and a very large number 
of levels (205) in the lower state, i.e. all bound levels 
(0 < v" < 48) and discretized scattering levels in a large 
energy range, with an energy up to 1100 cm -1 , described 
by physical or unphysical levels [17]. In this LbyL calcu- 
lation, the time-evolution of the total population in the 
excited state reproduces the one from the WP approach, 
in particular the low value of the population (~ 0.025) 
transferred at the end of the pulse. Furthermore, the 
time-dependence of the populations in the resonant level 
v' = 43 or in the level v' = 195 and also the variation of 
the total population in the bound 6 3 II levels, represented 
in the right panel of Fig. 5, reproduce perfectly the vari- 
ations calculated directly in the WP approach (Fig. 2). 
In the following, the set C is called the 'optimal' LbyL 
basis set. 

The wide energy range covered by the levels involved 
in the dynamics is not negligible compared to the fre- 
quency of the pulse Hojl ~ 7000 cm -1 . Therefore the 
validity of the PWA approximation is questionable. In- 
deed, for the pairs of levels E BjV >, E g . v » included in the 
basis set, the frequencies of the 'rotating' contributions 

[ "'" r — — u>i] are not always negligible compared to 

the frequencies of the neglected 'counter-rotating' contri- 
butions [ e '"' h g '" +u>l] (Appendix B 2). Further inves- 
tigation would be needed to check that the introduction 
of the counter-rotating terms does not change the main 
conclusions of the present analysis. 




FIG. 6: (Color online) LbyL description, optimal basis set C: 
variation of the population in the excited b 3 HV levels as a 
function of time (in ps). The laser pulse (vertical continuous 
and dashed blue [black] lines) satisfies the 7r-pulse condition. 
Total population of the excited bound levels v' — — 218 
(thick continuous black line). Non-adiabatic evolution for the 
population in the levels v' = 13 — 60 (thin continuous black 
line) resulting in a non vanishing final population transfer. 
Adiabatic evolution, following the variation of the pulse in- 
tensity, for the far-from-resonance levels v' = 167 — 218 (thin 
dashed red [black] line) located around the vertical transition 
v ~ 195 and the most populated at the maximum of the 
pulse or for the intermediate group of levels v' — 61 — 166 
(thin continuous orange [gray] line). 
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2. Two types of dynamics in the excited electronic state 

To go further in the analysis of the dynamics, we sep- 
arate the excited levels of the optimal set C into two 
different groups, according to the time-evolution of their 
individual population |6,/(£)| 2 . 

For levels < v' < 60 with a detuning varying in the 
range 2200 cm" 1 > 5 v > Vo » > —790 cm" 1 , the dynamics 
of population is very similar to that of the resonant level 
v' = 43. During the pulse the population |6„/(£)| 2 ex- 
hibits a small number of oscillations of a relatively small 
amplitude and some population remains in these levels 
after the pulse. The sum of the population in this group 
of levels grows almost monotonically during the pulse and 
reaches the final value ~ 0.02 (Fig. 6). 
As we move further off-resonant and consider bound lev- 
els 61 < v' < 218, we find that the evolution of the 
population is similar to that of the level v' — 195, i.e., 
traces time- variation of the pulse intensity ex [£(£)] 2 . The 
total population transferred at the maximum of the pulse 
t = tp is very high ~ 0.08 and it is larger than the popula- 
tion present in the group of levels close to the resonance. 
Yet no population remains after the end of the pulse. 
Thus it appears that two types of dynamics are observed 
for the levels of the excited electronic state. Levels with a 
not-too-large detuning remain populated after the laser 
pulse. Taking into account that the pulse is symmet- 
rical, Gaussian and unchirped, their evolution is neces- 
sarily non-adiabatic. Conversely, levels corresponding to 
highly-off-resonant excitation possess the largest popula- 
tion at the maximum of the pulse, but they do not retain 
their population after the pulse: such dynamics has thus 
a quasi-adiabatic character. Below we present a quali- 
tative description which emphasizes a relation between 
detuning and adiabaticity. 



B. Adiabaticity 

1. Introduction 

Adiabaticity of the evolution of a system is naturally ex- 
pressed in the basis of instantaneous eigenvectors of the 
Hamiltonian, the so-called adiabatic basic [20, 24] (see 
Appendix C). For a system with more than two lev- 
els, there is no general way to construct the instanta- 
neous adiabatic basis and thus no general expression of 
the adiabatic theorem [24]. In fact, the relationship be- 
tween adiabaticity, detuning, laser width and coupling 
strength can be perfectly illustrated in the case of a two- 
level system [v" , v'\ , where the instantaneous adiabatic 
levels can be explicitly constructed. The unperturbed 
vibrational levels \g) = \v") and |e) = \v'), define the 
diabatic basis (see Appendix C). The time-dependent 
wave function can be decomposed on the diabatic levels 
$!(t) = a(t)\g)+b(t)\e). We assume that only the level \g) 
is initially populated. The levels are coupled by a Gaus- 
sian pulse £of(t), with bandwidth Sui = 5.59 10" 4 au, as 



TABLE I: Different two-level systems [v" , v'] considered here. 
8 is the detuning of the pulse, (v'\v") is the overlap integral 
between the wave functions of the two levels and £l(tp) is the 
coupling at the maximum of pulse. The pulse has a duration 
tl — 0.12 ps and a bandwidth 8lo — 5.59 x 10" 4 au. 



system 


excitation 


8 [au] 


(v'\v") 


Q(tp) [au] 


a 


quasi 
resonant 


io-« 


0.03462 


2.10 x 10" 4 


b 


nearly 
resonant 


10" 4 


0.03462 


2.10 x 10" 4 


c 


out of 
resonance 


io-- 1 


0.03462 


2.10 x 10" 4 


d 


out of 
resonance 


2 x 10" d 


0.03462 


2.10 x 10 -4 


e 


quasi 
resonant 


io-* 


0.427364 


2.59 x 10" J 


/ 


far from 
resonance 


0.0198 


0.427364 


2.59 x 10 _d 



described in Section II B. In the RWA approximation, the 
time-dependent coupling is Q(t) = D x £$ (v'\v") f(t)/2. 
Below we study six different cases, labeled a to /; 
these differ by overlap integrals and dctunings (see 
Table I). For the overlap integral, we choose val- 
ues corresponding either to the resonant transition 
(v' Q = 43|uo" = 37) (systems a to d) or to the verti- 
cal transition (v' = 195|v " = 37) = 12.3 {v' \v ") (sys- 
tems e and /). The amplitude of the electric field So 
and the dipolc transition moment D are chosen such 
as, for (v'\v") = (v' \vo"), the n pulse condition, or 
Wl = VL£, is satisfied except for the cases e and /, 
where Of ~ 12.3 it. Therefore the maximum coupling is 
either smaller, Vl(tp) — 8uj /2.66 (cases a to d), or larger, 
fJ(tp) = 4.63 6u> (cases e, /), than the pulse bandwidth. 
We consider four classes of detunings. When the de- 
tuning satisfies \6\ <C 6u>, the systems are called 'quasi- 
resonant' (cases a and e). Systems where the detuning 
is | S | < 6ui are called 'nearly-resonant' (case b). Larger 
detunings such as |<5| > Su) correspond to 'off-resonant' 
excitation (cases c and d). 'Far- from- resonance' excita- 
tion, such as \6\ ^> Sui is represented by the case /, where 
the detuning S — 35.4 8ui is the one of the vertical tran- 
sition vq" — >• v' = 195. 



2. Two-level system: diabatic basis 

For the six considered cases, the time-variation of the 
population Pdiab{t) = |6(t)| 2 of the excited diabatic level 
|e) is calculated by solving the coupled system Eq. (C13). 
The results are presented in the first row of Fig. 7. For 
quasi-resonant systems, the time-dependence of the pop- 
ulation exhibits oscillations very similar to the Rabi os- 
cillations of a resonantly excited two-level system. In the 
case a, the population, initially in the ground level, is 
transferred continuously to the excited level during the 
pulse. In the quasi-resonant case e, the population oscil- 
lates more than 6 times from to 1 between the ground 
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and the excited levels, in agreement with the increase 
of the pulse area; at the end of the pulse 26.62% of the 
population remains in the excited level. For the nearly- 
resonant case b, almost all the population, up to 91.6%, 
is transferred monotonously to the excited level. 

For the off-resonant cases c and d and for the far-from- 
resonance case /, the time-evolution of the population 
in the excited level follows a Gaussian evolution similar 
to that of the pulse intensity ex \f(t)] 2 . After the end 
of the pulse, the total population returns to the initial 
level. When the maximum value of the coupling is small 
(for approximately Q(£p)/|<5| < 3/10), the perturbative 
limit is valid and Pdiabit) ~ [^{t)/S\ 2 ■ This occurs in the 
systems c, d and /, where Q(tp)/5 = 0.21, 0.10 and 0.13 
respectively; this can be also verified by comparing in 
Figs. 7c,d,f the first and fifth rows. The maximum value 
of the population in the excited diabatic level is small 
and varies proportionally to \(v' Q \vo")/5\ 2 . 



equation in the adiabatic basis (Eq. (C21)): 

Qa dl ab{t) = ^fi$(t)/ (E+(t) - E-(t)) 

1 5Q(t) 



2 [<5 2 + 4(f>(t)) 2 ] 3 / 2 



(21) 



In these equations, E±(t) (Eq. (C16)) denote the energies 
of the instantaneous adiabatic levels and 9 it) (Eq. (C15)) 
is the rotation angle occuring in the unitary matrix defin- 
ing the adiabatic levels. Here and hereafter, the dot in- 
dicates the time derivative. 

A strongly non-adiabatic instantaneous transfer corre- 
sponds to a high Qadiabit) value. From Eq. (21), one 
can deduce that non-adiabaticity occurs for a small de- 
tuning \5\, for a small coupling strength O(i), i.e. at the 
beginning and at the end of the pulse or for a pulse with 
low intensity, and also when the rotation angle 9(t) varies 
rapidly. 

The time-dependences of the rotation angle 6 '(£) /it and 
of the parameter Qadiabit) are shown respectively in the 
third and fourth rows of Fig. 7. 



3. Two-level system: adiabatic basis 



To resume the analysis of the adiabaticity of the pop- 
ulation transfer, we introduce now the adiabatic basis, 
made of the instantaneous eigenstates |\&±(i)) of the 
Hamiltonian (see Appendix C 4) . The time-dependent 
wave function is decomposed on the adiabatic levels 
*(£) = exp(-ift) [a(t)*_(t) + /3(t)*+(t)]. For a fully 
adiabatic process, the populations |a(i)| 2 and |/3(£)| 2 re- 
main constant during the pulse. In particular, if the sys- 
tem is initially in the adiabatic level |\&_(t = 0)) = \g), it 
remains in the instantaneous adiabatic level |\&_(£)) dur- 
ing the pulse and is in the level |^"-(i —¥ +oo)) after the 
end of the pulse. For a non-adiabatic process, the popu- 
lation of the adiabatic levels varies, and the stronger the 
non-adiabatic instantaneous population transfer is, the 
more rapid and important the changes of the instanta- 
neous adiabatic populations. In our study, the detuning 
S is fixed and only the coupling Q(t) varies with time. 
Therefore |v£_(£ — > +oo)) = \g). Simply, unless there is a 
significant non-adiabaticity, there is no population trans- 
fer. Conversely, a measure of the global non-adiabaticity 
of the process is related to the population transfer to the 
excited level after the pulse. 

The time-dependence of the population P a diab^) = 
|/3(£)| 2 in the adiabatic level |\&+(£)) is obtained by solv- 
ing the coupled system Eq. (C22), assuming that the pop- 
ulation is initially in the adiabatic level |\I/_ it = 0)) = \g) 
i\a{t — 0)| = 1). The results are drawn in the second row 
of Fig. 7. 

To characterize the adiabatic character of the instanta- 
neous population transfer, we introduce the parameter 
Qadiabit) deduced from the time-dependent Schrodingcr 



4- Two-level system: from quasi-resonant to 
far-Jrom-resonance excitation 

The described two-level model is useful for understanding 
the role of adiabaticity for both resonant to far-from- 
resonance excitations. 

For strictly resonant excitation, 6 = 0, and the rota- 
tion angle (Eq. (C17)) is equal to 6 it) — ^ at every 
time; then the initial wave function \g) corresponds to an 
equal mix of the two adiabatic levels \g) = l/\/2[|\IJ_(t = 
0)) — |\&+(i = 0))]. During the evolution, there is no non- 
adiabatic coupling (0 it) = 0) and no change in the pop- 
ulation of the adiabatic levels. The population oscillates 
between the ground and excited levels at the instanta- 
neous Rabi frequency Q(v',v" , t)/h. 
For quasi-resonant excitation, with a very small detuning 
\5\ <C £}(£p), as for example 5 = 10~ 8 au (cases a and 
e), the rotation angle 9(t) is almost equal to zero at the 
beginning and at the end of the pulse, when fi(t) <C \S\ 
(Fig. 7). Conversely, when Q(t) S> \6\, the rotation angle 
remains constant and equal to 9 = ? . For t ~ t na± , with 
2\n(t na± )\ = \6\ (or 9it na± ) = §), the rotation angle 
changes rapidly. Two sets of nearly adiabatic levels can 
be introduced, \&±(£), valid at the beginning t < t na _ 
or at the end of the pulse t > t na+ and v£±(t), valid 
during the pulse t na _ < t < t na+ . During these three 
time intervals, the evolution is completely adiabatic. For 
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adiabatic couplings Qadiabit na± ) ■- 2STT V LLL ^~\s\ 
huge, of the order of 3 x 10 4 in the quoted examples. 
Therefore, strong instantaneous population transfer be- 
tween the adiabatic levels occurs only around t na± ■ The 
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FIG. 7: (Color online) Two- level system: analysis, as a function of time (in ps), of the dynamics for the different cases a to f 
described in Table I. The maximum and the duration of the pulse are indicated by the vertical continuous and dashed blue 
[black] lines. First row: population Pdiabit) in the diabatic excited level [e), the initial population lying in the lower diabatic 
level \g). Second row: variation with time of the population P^ diab (t) in the adiabatic level ^/+(t), the initial population lying 
in the adiabatic level ^_(t = 0) = \g). Third row: rotation angle 9{t)/ir (Eq. (C17)) defining the instantaneous adiabatic 
levels $±(i). Fourth row: parameter Q a diab(t) (Eq. (21)) characterizing the adiabaticy of the instantaneous transfer. Fifth 
row: pulse parameter Q(t)/5. The quantities have been multiplied by the factor indicated in the upper left corner. 



value \b na ± — tp\ ~ 0.34 ps is much larger than the pulse 
duration tl — 0.12 ps. The population transfer occurs 
thus in the wings of the Gaussian pulse, at the turn-on 
and turn-off of the pulse, when the laser intensity is al- 
most negligible. For t < t na _ , the population remains in 
the lowest adiabatic level, described by the wave function 
\£?(t) = \g)exp[—iE g i\. If one sets E g = 0, then there 
is no change in the phase of this wave function. For 



t na - < t < t na+ , 9{t) = ~ and the adiabatic levels "J/ 



±(«) 



correspond at each time to an equal mix of both diabatic 
levels with a phase varying with time (Eq. (C25)). These 



adiabatic levels evolve as 



m 



it)) 



V2 



[T\g) + \e)\ exp[- 



E±(t')dt'] 



For tna- < t < tna+i the wave function \I/(£) can be 
decomposed on these states, with amplitudes a p (t) and 
(3 p (t). The absolute values of these amplitudes remain 
constant (see row 2 in Fig. 7a, e). These populations 
can be estimated in the sudden approximation [24], by 
projecting the adiabatic wavefunction \ 1 $' _(t na _)), valid 
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just before t na _, on the adiabatic functions ^>^.{t na _), 
valid just after t na _. In this way, one obtains |a p (£)| = 
|/3 P (£)| = -t=, just like for the strictly resonant excita- 
tion. The Rabi oscillations occurring at quasi resonance 
(Fig. 7e) in the population of the excited diabatic lev- 
els | < ^>(t)|e)| 2 , for t na _ < t < t na+ , result from a 
beating effect in the coherent superposition of the *&±{t) 
adiabatic levels. At t = t na+ , a strong non-adiabatic 
coupling occurs again during a very short time. The 
sudden approximation allows one again to obtain the 
value of the final population of the diabatic levels \g) 
and |e) after the end of the pulse, in terms of the quan- 
tity 9' = f t ™ + [E + (t) - E_{t)]dt, which, in the limit 
S — > 0, is equal to the pulse area 9(+oo). 
With increasing S, \t na± -t P \ and Qadiab(t na± ) decrease, 
being equal in the nearly- resonant system b to \t na± — 
tp\ ~ 0.15 ps and Q a diab(t na± ) ~ 1.2. Non-adiabatic 
transfer of population from the adiabatic level <3/_ (t) to 
the adiabatic level ^+(t) occurs in two steps around t na± , 
but with a less-pronounced sudden character. 
The maximum value of 0(tp) decreases when \S\ in- 
creases, and for a detuning such as \S\ > 2Q,{tp) or 
0(tp) < j, as in cases c, d and /, the times t na± do not 
exist. The maxima of the parameter Q a diab(t) become 
smaller and appear during the pulse at t± (\t± — tp\ < 



tl ) , with t± = tp±- 



-Mi- 



4n(t P ' 



1. As a result, the 



transfer becomes more adiabatic, with a very low popu- 
lation Padi a b(t) transferred to the adiabatic level '!'+(£) 
at t ~ £_ . In addition, this population transferred to the 
upper adiabatic level ^S+{t) returns back to the lower 
adiabatic level ^~{t) at t ~ t+. Almost no population 
remains in the excited level after the pulse. 
For a sufficiently high \6\ value, the evolution of the 
population transfer becomes completely reversible and 
the population in the excited adiabatic level is such as 
\f3(tp — St)\ = \(l{t P + 6t)\, as observed in systems d 
and /. The population of the adiabatic excited level can 
be calculated in the perturbative approximation, leading 

[Qadiab(t)] 2 (compare 
For a very large 



■.,.., 



t0 "odiabW ~ \-2[E+(t)-E-(t)]l 

rows 2 and 4 in Fig. 7 for system d). 

detuning, \5\ ^> fi(tp), the population of the adiabatic 

excited level is maximum at t = t± with 
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(22) 



For off- resonant excitation, |<5| > Slu, the maximum pop- 
ulation in the adiabatic level is very small {PgtLt.it ±) <C 
1), decreasing with \6\ more rapidly than the maximum 
population in the diabatic levels, equal, in the perturba- 
tive limit, to: 



characterizing the excitation of a two-level system by 
a Gaussian pulse fully determine the dynamics. The 
nearly-resonant or off-resonant character of the process 
depends on the ratio S/Sco. For a nearly-resonant ex- 
citation of a level |t>') lying within the pulse band- 
width {\S v i jVo " | < 5lu) and in the weak field limit 
fi(tp) <C 5ui, the population transferred to this level 



is \b v ,{t -»• ^)| 2 = 8^1n2[^ij [^j . Nearly- 
resonant excitation acquires a strong-field character as 
soon as Ct{tp) ^> Sco; then a highly non-adiabatic trans- 
fer occurs in the wings of the pulse, at t na ± . The final 
population transfer strongly depends on the value of the 
phase difference accumulated in the adiabatic wavefunc- 
tions ^±{t) during the time interval [t na -, t na +}- This 
phase difference is at the origin of the Rabi oscillations 
observed in the population of the diabatic levels \g) and 
|e>. 

For off-resonant excitation, the evolution can be de- 
scribed in the perturbative limit if Cl{tp) <C 5; it re- 
sults in a completely adiabatic dynamics, with no final 
population transfer, the population of the diabatic levels 
following the variation [/(t)] 2 of the pulse intensity. 
Let us emphasize that the conclusions reported above are 
valid for a pulse with sufficiently slow time-dependences 
in both the electric field envelop and the instantaneous 
frequency In particular, they are not valid for a spec- 
trally cut Gaussian pulse [25]. In this case, rapid vari- 
ations of the instantaneous frequency around the pulse 
maximum are responsible for a nonadiabatic character 
of the off-resonant excitation and the off-resonant levels 
remain populated after the pulse [26]. 



5. Excitation of molecular wavepackets 

Some information on the adiabaticity of the population 
transfer in the optimal multi-level system can be ob- 
tained by considering directly the excitation of molecular 
wavepackets in the lower and excited electronic states and 
reducing the WP description to a two-level problem. In- 
deed, as we are considering pulses with duration (~ 100 
fs) much larger than the vibrational period of the consid- 
ered levels, it is possible to ignore the vibrational motion 
during the laser excitation, that is to study the excitation 
process in the impulsive approximation [27]. 
In the time-dependent Hamiltonian describing the laser 
excitation of a diatomic molecule in the WP approach, 
given in Eq. (B2), 



H = 



T + V g {R) + 
\W eg {t) 



2 



W ge {t) 



T+V e {R)- 



2 



,(24) 



Pdiab{t} 



n{t P ) 



(23) 



we neglect the kinetic energy T. Introducing the differ- 
ence between the two dressed potentials 



To summarize, the three parameters 5, Slu and Q{tp) 



A L {R) = V e {R)-V g {R)-ftujL 
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FIG. 8: (Color online) Criteria of adiabaticity for the excita- 
tion of a molecular wavepacket by a gaussian pulse: variation 
of Q a diab(R>t) (Sec. IV B 5) as a function of the internuclear 
distance R (in au) and of the time t (in ps). The condi- 
tion for an adiabatic evolution is broken near R c = 10.5 au 
and for times corresponding to the beginning and the end of 
the Gaussian pulse. Non-adiabatic effects are more impor- 
tant for low-intensity pulse W^/120 (upper panel) than for 
high-intensity pulse WZ (lower panel). 



and ignoring the mean potential 

2V(R) = V e (R) + V g (R), 

which introduces only an independent phase factor, we 
can consider, at each internuclear distance R, the two- 
level Hamiltonian in the diabatic representation [8] , 



, 2 



A L (R) 

,(*) 



\w ge (t) 

f±A 



(R) 



(25) 



and analyze the adiabaticity of the excitation process, by 
calculating the R— and ^-dependent function Q adiab (R, t) 
similar to the function Q a diab(t) defined in Eq. (21). 

Figure 8 shows Q a diab{Ri *) as a function of R and t. For 
a pulse with a carrier frequency ul resonant with the 



transition a 3 S + 



vo 



37 -> b 3 nv' n 



43, the dressed po- 
tentials cross each other (Al(R c ) — 0) at R c = 10.53 
au. The adiabatic condition is broken at the beginning 
and at the end of the Gaussian pulse t = tp — tl = 0.48 
ps and t — tp + tl = 0.72 ps (small laser intensity). It 
is also broken at internuclear distance close to R c (small 
detuning). This determines both the times where popu- 
lation can be transferred from the ground electronic state 
to the excited electronic state and the spatial location of 
the transferred population. For the studied unchirped 
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FIG. 9: (Color online) LbyL calculation: influence on the 
dynamics of the far-from-resonance b 3 Hv' levels, when only 
the level vo" = 37 is introduced in the lower state. Variation 
with time (in ps) of the populations of the vibrational levels 
6 3 I1 v' with v' = 43 (thick continuous orange [gray] line), v' = 
42 (thin continuous orange [gray] line), v' — 44 (thin dashed 
orange [gray] line), v — 195 (thick black line), v — 194 (thin 
black line), v' = 196 (thin dot-dashed black line) and of the 
initially populated vibrational level a E + vo" = 37 (dashed 
red [black] line). Some populations have been multiplyied by 
the factor indicated in the figure. The pulse maximum at tp 
(vertical continuous blue [black] line) and duration indicated 
by the vertical dashed blue [black] lines satisfies the 7r-pulse 
condition. Different basis sets [«" , v'\ are used. The difference 
between lefthand and righthand panels of the same line is 
the absence or presence of far-from-resonance levels. Upper 
line, a) [vo" = 37, v' = 43] and b) [v " = 37, v' = 43 and 
v' = 194 - 196]; middle line, c) \v " = 37, v' = 42 - 44] and 
d) [vo" = 37, v' =43- 44 and v' = 194 - 196]; lower line, e) 
[vo" = 37, v' = 13 - 60] and f) [v " = 37, v' = - 218]. 



pulse, population transfer occurs around R c and in the 
wings of the pulse. 



Influence on the dynamics of the 
far-from-resonance levels 



We pay now a particular attention to the contribution 
to the dynamics coming from far-from-resonance levels. 
We first analyze the dynamics of excitation in a multi- 
level system including both nearly- resonant and far-from- 
resonance excited levels bv' while keeping only the single 
level avo" in the lower electronic state. Then we intro- 
duce all the lower [v" — — 204] levels of the optimal 
set to obtain a complete view of the modifications of the 
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end of the pulse. In fact, for the two ^o" an d v' lev- 
els, which are degenerate in the diabatic representation 
(Eq. (C3) with A e v , =43v , = Eq. (11)), the maximum 
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FIG. 10: (Color online) LbyL calculations: influence on the 
dynamics of far-from-resonance 6 3 II v' excited levels, when 
all the levels a 3 £ + v" of the set C are introduced in the 
lower state. Basis sets [v" = — 204, v' — 13 — 60] (orange 
[gray] lines) and [v" = - 204, v' = - 218] (red [black] 
lines). The pulse (vertical continuous and dashed blue [dark] 
lines) satisfies the 7r-pulse condition. Left column: variation 
with time (in ps) of the population of the 6 3 II v' levels (in 
percentage); panel a): total population in the levels v' — 
13 — 60; panel b) population in the resonant level v' — 43. 
Right column: final distribution of population (in percentage) 
in the excited levels v' (panel c)) and lower levels v" (panel 
d)) with a maximum, for vo" , equal to 90.9 and 76.0 for the 
first and second basis sets respectively. 



dynamics of close to the resonance levels induced by far- 
from-resonance levels. 



1. A single level in the lower electronic state 

Starting with a single vq" = 37 level in the ground elec- 
tronic state, we progressively grow the basis set in the 
excited state, by adding either nearly-resonant levels, i.e. 
close to the oblique transition, or far-from-resonance lev- 
els, i.e. close to the vertical transition. 
In the upper row (Figs. 9a, b), we compare the two- 
level system that consists only of the two resonant levels 
[vo" = 37,v' = 43] to the small 5-level system [vq" — 
37, v' = 43, 194 — 196] containing the two resonant levels 
and the three far-from-resonance levels v' = 194 — 196. 
As expected, for the resonant two-level system excited 
by a 7r-pulse, a total exchange of population is observed. 
For the basis set b), a very low population is trans- 
ferred during the pulse to the three additional levels, 
less than 1%. Nevertheless the presence of these lev- 
els modifies completely the dynamics of the population 
in the v' — 43 resonant level: only 29% of the popu- 
lation is transferred in this level at t — tp, instead of 
50% in the two-level system, and, in addition, the pop- 
ulation disappears almost completely (0.012%) at the 



coupling strength, Cl v 



0=43, i>o' 



= 2.10x 10" 4 au, is of the 



same order of magnitude as the second order contribution 
[^v=i95,v ") 2 / Ay, =195 y o = 3.4 x 10~ 4 au corresponding 
to the vertical transition v' — 195. This modifies signif- 
icantly the energies of the instantaneous adiabatic levels 
connected to the resonant vq" = 37 and v' = 43 dia- 
batic levels and therefore the phase difference accumu- 
lated between t na _ and t na+ in the adiabatic wavefunc- 
tions fy± (Eq. (22)) or, during the pulse, the beating 
between the probability amplitudes of the resonant lev- 
els (Sec. IV B 4). This explains qualitatively the strong 
changes in the dynamics of the excitation process. 
In the middle row (Figs. 9c, d), the two nearly-resonant 
levels v' = 42 and v' = 44 are added to each above de- 
scribed basis set. For these levels, the overlap integrals 
with the wo" level (0.029 and 0.031) are nearly equal to 
the overlap integral (0.034) of the resonant transition, 
and the detunings are small ±0.0002 au. In the system 
[vq" , v' — 42 — 44] introducing only the nearly-resonant 
levels, the population transferred to the excited state is 
shared between the three excited levels with a large to- 
tal transfer (82%). When the three levels close to the 
vertical transition v' = 194 — 196 are added, there is a 
low transfer (3.3%) to the resonant v' = 43 level but a 
larger transfer (12.5 and 15.6%) to the nearly-resonant 
levels v' = 42 and v' = 44. Here also the introduction of 
the far-from-resonance levels modifies the excitation dy- 
namics of the nearly-resonant levels. In particular, there 
is a strong decrease in the population transferred to the 
excited electronic state at the end of the pulse, 31.4% in- 
stead of 82.3%, and an important change in the branching 
ratios in the population of the nearly-resonant levels. 
In the lower row (Figs. 9e,f), larger basis sets are intro- 
duced in the excited state. The set e) [«o" = 37, v' — 
13 — 60] consists on all the excited levels remaining pop- 
ulated after the pulse (Fig. 6) in the WP calculation. 
The dynamics of excitation of the nearly-resonant lev- 
els v' — 42 — 44 is qualitatively the same as in set c), 
with a total population transfer equal to 78.2%, but with 
a change in the branching ratios. This shows that, in 
this group of v' = 13 — 60 levels, the nonadiabatic dy- 
namics (Sec. IV A 2) is dominated by the three nearly- 
resonant levels. For the set f) [v " — 37, v' = — 218] 
introducing a still larger basis in the excited state, the 
dynamics of the excitation of the nearly-resonant lev- 
els differs from that one observed for the set d) [«o" = 
37, v' = 42 - 44, 194 - 196]. The final population in 
the nearly-resonant levels v' — 42 — 44 decreases from 
the value 31.4% to the value 13.5%, showing that far- 
from-resonance levels other than the v' — 194 — 196 ones 
contribute to the dynamics. For the two basis sets e) 
and f), almost all the population is transferred to the 
excited electronic state, 93.1% and 98.24% respectively 
(see the low value of the final population of the initial 
level v$" = 37), this population being mainly distributed 
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in the nearly- resonant levels v' ~ 42 — 44 for set e) , but 
in lower levels v' ~ 36 — 40 for set f ) . 



2. Several levels in the lower electronic state 

In this section we incorporate all the levels v" = — 204 of 
the lower electronic state included in the optimal basis set 
(Sec. IV A 1) and we analyze how the far-from-resonancc 
levels modify the excitation dynamics. We consider the 
basis sets [v n = - 204, v' = 13 - 60] and [v" = - 
204, v' = — 218]. The first set encompasses only the 
excited levels which remain populated after the pulse in 
the WP treatment whereas the second one is the optimal 
set. In Fig. 10 we show the computed time-evolution of 
the total populations of the levels v' = 13 — 60 (Fig. 10a) 
and of the resonant level v' = 43 (Fig. 10b). We display 
the final distributions of population in the vibrational 
levels of the excited (Fig. 10c) and the ground (Fig. lOd) 
electronic states. 

At the beginning of the pulse, for t < 0.53 ps, when 
the pulse intensity increases, weak population recycling 
has yet occured and the contribution to the dynamics of 
far-from-resonance levels is not very important. When 
the pulse reaches its maximum intensity at t = tp = 0.6 
ps, the population exchange between the lower and ex- 
cited electronic states becomes more important and no- 
ticeable changes become observable in the Rabi oscilla- 
tions occurring either in the population of the resonant 



to what is observed in Fig. 5 using the optimal set. Fur- 
thermore there is a strong transfer of population to the 
excited electronic state (98.2%), substantially different 
from the weak transfer (1.9%) obtained with the optimal 
set. In this section, we analyze in more detail the specific 
role of the levels of the ground electronic state, especially 
those which are quasi-degenerate with the initially pop- 
ulated one. In the energy range close to this initial level, 
the spacing between consecutive bound levels is indeed 
much smaller than the laser bandwidth Scu = 120 cm -1 . 
All bound levels with v" > 34 and all continuum states 
with an energy up to 115 cm -1 (with the chosen grid, dis- 
cretized scattering levels up to v" = 134) are such that 
|A* W |<<^. 

To analyze the excitation dynamics from a group of 
quasi-degenerate levels, we consider in Appendix D a sim- 
ple model describing the excitation of a single sublevel 
from an iV-fold degenerate level, which admits an ana- 
lytical solution. The comparison with a non-degenerate 
two-level system shows that, in the high field regime, 
the population transfer is divided by N, whereas the re- 
maining population is equally distributed among the N 
sublcvcls of the degenerate lower level. When the num- 
ber N increases, the transfer of population toward the 
excited level decreases: there is a blockade of the excita- 
tion induced by the degeneracy of the lower level, with 
no transfer at all for N — > oo. 

As the first example we consider two basis sets, [v" = 
36 - 38, v 1 = 43] and [v n = 36 - 39, v' = 43], and the 



level Vq — 43 or in the total population of the levels 7r-pulse resonant with the transition v" = 37 



v = 13 — 60. When far-from-resonance levels are intro- 
duced in the basis set, the total population transferred 
during the pulse to the levels v' = 13 — 60 is smaller; 
nevertheless there is no significant change in the total 
population transferred to the excited state (1.9% instead 
of 1.8 %)(Fig. 10a). Concerning the Rabi oscillations of 
the population of the resonant level v' = 43 (Fig. 10b), 
some modifications occur for t > tp and half of the pop- 
ulation remaining in this level is transferred back to the 
ground state (0.04% instead of 0.11%). We remind the 
reader that a similar decrease in the population of the 
resonant level induced by including far-from-resonancc 
levels has already been observed in the simple LbyL mod- 
els discussed in Fig. 9). Concerning the final distribution 
in the excited state population(Fig. 10c), the population 
is, on average, shifted toward slightly higher w'-values. In 
the ground state (Fig. lOd), the population is spread over 
a larger energy domain, with a smaller population trans- 
ferred back to the initial level uo" = 37 (76.0% instead of 
90.9%). 



D. Blockade of the excitation due to a 
quasi-degenerate level in the lower electronic state 

When the basis set [vo" = 37, v' = — 218] is used 
(Figs. 9e,f), the resonant level v' — 43 population does 
not exhibit a large number of Rabi oscillations contrary 



43 (Figs. lla,b). For the set a) where N=3, the final 
population transferred to the v' level is equal to 0.0556. 
The populations of the initial t>o" level and of the other 
two ground levels are respectively equal to 0.1313 and 
0.4065. For the set b) with TV = 4, no population remains 
after the pulse in the excited level v' , the effective pulse 
area being equal to ^j^-tt — n. Simultaneously, the four 
ground state levels are equally populated (population j). 
At the maximum of the pulse, when -^6(ip) = \ , the 
population in the excited level is equal to -tt = \. In 
the initially populated level t>o" it is equal to yg and 
in the other three ground levels to A. The evolutions of 
the populations in the quasi-degenerate case agree almost 
perfectly with the iV-fold degenerate model with TV = 3 
and N = 4. 

More realistic results, obtained by using the LbyL ap- 
proach with the basis set [v" — 32 — 42, v' = 43], are 
also presented (Figs. llc,d,e), together with comparable 
results of the analytical model. For the ground levels 
|A^„ v „ | < 0.000027 au and the overlap integrals vary in 
the range 0.018-0.040. For this system, which includes 
a 11-fold quasi-degenerate lower level, the dynamics of 
excitation is similar to that of a degenerate level with 
N = 11. The effective area for the pulse is y/Nn = 3.37T, 
and, during the pulse, the population of the resonantly 
excited level oscillates between and -^ = 0.09. Simulta- 
neously, the population is redistributed among the quasi- 
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FIG. 11: (Color online) Influence of a quasi-degenerate level in 
the ground electronic state. Evolution of the populations as a 
function of time (in ps), under a pulse (continuous and dashed 
blue vertical lines) satisfying the 7r-pulse condition W = W£. 
The results obtained using the LbyL method are compared to 
the analytical formulas (Eq. (D5)) describing the excitation 
from a JV-fold degenerate level (thin continuous red [black] 
lines), where the initial population is in a single level (AJ of 
the manifold, the other N — 1 sub-levels (A.^) being unpopu- 
lated; the excited level (Bj) is resonantly excited. Panel a) 
N = 3: basis set [v" = 36 - 38, v' = 43]. Panel b) N = 4: 
basis set [v" = 36 - 39, v' = 43]. Panels c), d), e) N = 11: 
basis set [«" = 32 — 42, v' — 43]. Initially populated level 
a 3 E + vo" = 37 (double-dot dashed black line) and resonantly 
excited level b 3 Ilv' = 43 (thick continuous clear-gray line]. 
Non resonant a 3 £ + levels described by the same formula in 
the analytical model: v" = 36 (medium-thick dot-dashed or- 
ange [gray] line); v" — 38 (medium-thick dashed orange [gray] 
line); v" = 39 (medium-thick continuous orange [gray] line). 



degenerate levels of the ground electronic state. 
To summarize, in the strong-field excitation from a level 
close to the dissociation threshold, the high density of 
levels in the initial state is at the origin of a blockade of 
the excitation process. Simultaneously the increase of the 
effective Rabi frequency explains the oscillations occuring 
during the pulse in the population of the quasi-resonant 
levels. Let us remark that this phenomenon is similar to 
the ionization suppression occuring in the Rydberg atom 
ionization by an intense laser pulse. When the N initial 
discrete levels are exactly degenerate, only I/TV of the 
initial population ionizes in a time divided by the factor 
N [28]. 

Finally we remark on the influence of the far-from- 
resonance excited levels on the excitation from quasi- 
degenerate lower levels. As expected, the basis set 
[v" — 0—204, v' = 43] yields a blockade of the excitation: 
the final population of the initial level t>o" = 37 is large, 
amounting to 90.6% and, simultaneously, a very weak 



population, equal to 1.0%, is transferred to the resonant 
level v' — 43. When adding some far- from- resonance 
excited levels, in the basis set [v" = — 204, v' — 
43, 194— 196], one observes simultaneously a blockade of 
the excitation and an important redistribution of popula- 
tion within the ground state: the population transferred 
to the resonant level v' = 43 is almost negligible (1%), 
but 78.5% of the population is redistributed among the 
levels 20 < v" < 47 and 20.5% in the continuum, mainly 
in scattering levels with an energy smaller than 20 cm . 
Here also, the contribution of the far-from-resonance lev- 
els v' = 194— 196 is crucial. These levels are only weakly 
populated during the pulse, but their population is recy- 
cled back to a large number of vibrational levels of the 
ground electronic state. 



V. DISCUSSION AND PERSPECTIVE; TRAIN 
OF PULSES 

In this paper, we have explored the possibility of enhanc- 
ing the rate of formation of stable RbCs molecules in the 
absolute ground level Rb(5s)Cs(6s) X^+v" = 0. More 
precisely, we have analyzed the excitation by a single 
unchirped Gaussian pulse of molecules already formed 
in weakly bound levels of the Rb(5s)Cs(6s) a 3 E + state, 
after photoassociation of ultracold Rb and Cs atoms fol- 
lowed by spontaneous radiative decay. When the final 
level of the excitation is a spin-mixed level, for example 
a level of symmetry + or 1, it will be possible to trans- 
fer optically, in a second step, the population from this 
electronic excited level to the absolute ground level. Re- 
stricting the description to uncoupled electronic states in 
the Hund's case a coupling scheme, we have investigated 
the possibilities offered by the presently widely developed 
femtosecond laser sources, to transfer efficiently the pop- 
ulation from the excited a 3 £ + t>o" = 37 level toward the 
b 3 Uv' = 43 level. 

The dynamics of the photoexcitation process is mod- 
eled with the WP method, calculating the time-evolution 
of the wavepackets propagating along the electronic 
states coupled by the laser pulse. This employed non- 
perturbative method allows us to analyze dynamics all 
the way from the low- to the high-intensity regime, the 
latter being the 7r-pulse for the resonantly-driven vq" = 
37 — > v' = 43 transition. We have also developed a LbyL 
description, and employed a variety of subsets of vibra- 
tional levels in numerical simulations. The comparison 
between various restricted LbyL calculations and with 
the full-scale WP results allowed us to qualitatively un- 
derstand the complex evolution of the population of the 
multitude of vibrational levels and to identify the specific 
influence of these levels on the dynamics. 
In the perturbative limit, only the quasi- resonant levels 
lying within the bandwidth of the intensity distribution 
are excited. For a not too short duration of the pulse T£ 
corresponding to a bandwidth Suj smaller than the energy 
spacing of the vibrational levels in the excited electronic 
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state, the excitation process is selective. Conversely, the 
efficiency of the population transfer is very low. 

With increasing intensity, more levels are excited, but 
the excitation rates remain very small and the dynam- 
ics becomes more adiabatic. We have shown that, in 
the strong-field regime, i.e. for the 7r-pulse, far-from- 
resonance excited levels with high Franck-Condon fac- 
tors are populated during the pulse. At the maximum of 
the pulse intensity, the population of the levels excited 
through a vertical transition at the outer turning points 
of their wave functions is much larger than the popu- 
lation of the quasi- resonant levels. Because of the adi- 
abatic character of the excitation of far-from-resonance 
levels, the time-dependence of their population follows 
the smooth evolution of the pulse intensity, notwith- 
standing the high value of the Rabi frequency, and after 
the pulse no population is retained in these levels. Nev- 
ertheless these levels contribute to the excitation pro- 
cess, in particular they give rise to a large population 
recycling and to an important population redistribution 
in the ground electronic state. The individual contribu- 
tions to the amplitude of population of a particular level 
arising from the other levels are very intricate, making 
the analysis and therefore the control of the excitation 
process with unchirped pulses in the high-field regime 
very difficult. It is worth noticing that quantum con- 
trol of molecular wavepackets could still be possible in 
the strong- field regime [29], but in very specific condi- 
tions: for 'vertical' transitions with high Franck-Condon 
factors, not involving quasi-degenerate lower levels, and 
imposing special shape requirements to the standing edge 
of the pulse [30] , in order to avoid population redistribu- 
tion in the lower state. 

Furthermore, for the system under study, the initial level 
avo" — 37 lies very close to the dissociation thresh- 
old, where the density of levels is very high. We have 
shown that this situation is at the origin, in the high-field 
regime, of the excitation blockade. A simple model de- 
scribing the resonant excitation of a single level from an 
iV-fold degenerate level is developed, which reproduces 
comparable LbyL calculations well. In the strong-field 
regime, the excited level population is governed by a 
pulse area larger by a factor of y/N than the real area 
of the pulse and its amplitude is divided by N. This 
explains the oscillations with a low amplitude observed, 
during the pulse, in the population of the resonantly- 
excited level. Due to this blockade phenomenon, high 
intensity and large bandwidth pulses are poorly suited 
for gaining high excitation rates. 

In order to increase the population transfer to the reso- 
nant level v' — 43 while conserving the selectivity pro- 
vided by the low field regime, we have explored the possi- 
bility of using a weak intensity train of ultrashort coher- 
ent pulses. We recapitulate main features of such trains 
in the time and frequency domains in Appendix E. Pre- 
viously, coherent excitation of a two-level system by a 
train of short pulses has been described analytically [31]. 
Transient coherent accumulation for two-photon absorp- 



0.015 




01234567 12345678 



© 


: c) ' ' Jf"~- 




-d) s- 


- 


'■3 0.8 




_ / 


_ 


« 


* 


j 




ft 0-6 


/ 


1 


- 


© 


M 


j 


. 


ft 0.4 


/ ~ 


- j 


- 


0.2 


/ 


- J 


- 




Vi, , , i , , , i , , , i , , , - 


y i , , , i 


i . i 1 i i i 



0.015 

0.01 

0.005 


1 

0.8 

0.6 

0.4 

0.2 



20 40 60 80 20 40 60 80 

time [ ps ] 

FIG. 12: (Color online) Panels a), b), c): Excitation of the 
resonant transition a 3 £ + vo" — 37 — ► b 3 Tlv' = 43 in the 
low field regime by successive Gaussian femtosecond pulses 
with duration n, — 0.12 ps, repetition time T rep — 0.8 ps and 
maximum coupling strength W£/120. Variation with time (in 
ps) of the population of the excited electronic state 6 3 II (thick 
black line), of the resonant level b 3 Hv' = 43 (thin red [black] 
line) and of the quasi-resonant levels v' = 42 (thin orange 
[gray] line) and v' = 44 (thin black dashed line) levels. Panel 
a): train of TV = 10 coherent pulses in the LbyL approach with 
the basis set [v " — 37, v' = 42- 44]. Panel b): train of Af = 
10 pulses in the WP approach, where the lower and excited 
wavepackets obtained in the middle of the interval between 
the pulses Pi and P*+i are taken as initial wavepackets for 
the pulse Pi+i. Panel c): same as for panel a), but for a 
train of TV = 120 pulses, corresponding to an effective ir- 
pulse. Panel d): single picosecond 7r-pulse pulse with duration 
120tl = 14 .4 ps and maximum coupling strength Wl /120 in 
the LbyL approach with the same basis set. 



tion via an intermediate level has been demonstrated in 
atomic Rb [11]. Application to efficient selective vibra- 
tional population transfer between electronic states of a 
diatomic molecule has been discussed by Araujo [32]. 
In a preliminary study, we analyzed the dynamics of ex- 
citation by a train of pulses in the perturbative regime. 
Each individual pulse has the Gaussian shape of dura- 
tion T£ = 0.12 ps and of maximum coupling strength 
Wl/120. The repetition time is T rep = 0.8 ps, slightly 
smaller than the vibrational period T mb , in the excited 

r e,v 

state (Sec. Ill A) and with a vanishing pulse-to-pulse 
carrier-envelope-offset phase shift A<fi ce . We performed 
calculations both in the WP and the LbyL approaches. 
In the WP calculations, the final state at the end of each 
individual pulse is taken as the initial condition for the 
following pulse. The total population transferred to the 
excited electronic state b 3 Tl and its distribution among 
the different v' levels during Af = 15 pulses are shown in 
Fig. 12 (upper right panel). We also carried out LbyL cal- 
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culations using the basis set [vo" — 37, v' = 42 — 44]. The 
computed time- variation of the population in the excited 
levels is presented in Fig. 12 (left panels). The LbyL cal- 
culations reproduce perfectly the WP results, the small 
LbyL basis set clearly being sufficient for low-intensity 
pulses. 

Pulse after pulse, there is an accumulation of the total 
population transferred to the excited electronic state & 3 IL 
When the number of pulses, J\f, increases, the distribu- 
tion of population among the excited levels bv' becomes 
more selective, with an accumulation of population in 
the resonant level v' . A given pulse transfers to the level 
bv' an probability amplitude which interferes with the al- 
ready present probability amplitude, transferred by the 
previous pulses. The nature of the interferences depends 
on the phase exp[iT rep d v > tVQ »] involving the detuning of 
the considered bv' level [32]. For the resonant level bv' 
the interferences are constructive and the population in- 
creases with the growing number of pulses. For other 
vibrational levels, due to the mismatch in this phase, 
the population will oscillate with A/", without experienc- 
ing accumulation. This increase of the selectivity of the 
excitation with the number of pulses is a signature of 
the comblike structure of the energy spectrum £(ui) (Eq. 
(E2)) of the pulse train. The frequency spectrum consists 
of equally spaced "teeth" , with a spacing proportional to 
the repetition frequency T^l, an intensity proportional 
to A/" 2 and a width narrowing as Af~ 1 . 
When the number of pulses becomes equal to 120, the to- 
tal population initially in the vq v — 37 level is transferred 
to the resonant excited level v' Q =43. Let us remark that 
the total pulse area for the train with 120 pulses, each 
with a pulse area equal to 7r/120, amounts to tx. This 
train of pulses is thus equivalent to a single 7r pulse in 
the perturbative regime. If the number of pulses contin- 
ues to increase, the population cycles back to the initial 
level. The total duration of the pulse train amounts to 96 
ps, much smaller than the lifetime of the resonant level, 
which is smaller than 30 ns [17]. 

In conclusion, to succeed in controlling the dynamics 
of photoexcitation with unchirped femtosecond lasers, it 
seems necessary to employ low-intensity pulses. The ra- 
diative lifetime of the excited level can be disregarded 
if the total duration T to t — NT rep of the pulse train is 
sufficiently small. Consequently, the repetition rate has 
to be as large as possible. We notice the technological 
developments aimed at increasing the repetition rates in 
the train of femtosecond or a few picosecond pulses us- 
ing acousto-optic devices are in progress. Repetition fre- 
quencies up to 50 GHz were obtained with electro-optic 
phase- modulator shaping of a picosecond laser [33] , much 
faster than those of the order of 100 MHz obtained from 
a Kerr lens mode-locked femtosecond Ti-sapphire laser 
[11] or of 100 kHz for a regenerative amplifier seeded by 
a Mira oscillator [34] . 

For completeness, we mention here that another way 
to obtain high transfer rate with high selectivity re- 
lies on the use of a single pulse in the picosecond do- 



main with a sufficiently narrow bandwidth. To illus- 
trate this point, we computed excitation by a single 
Gaussian pulse, resonant with the transition «o" — * v 'oi 
with maximum coupling W' L = Wl/120 and duration 
t' l = 120 x tl = 14.4 ps (Fig. 12, lower right panel). 
Calculations were done in the LbyL method using the 
basis set [vo" = 37, v' = 41 - 44]. The bandwidth of this 
pulse, equal to Suj/120 = 1 cm -1 , is sufficiently narrow 
to include only the single v' level within its bandwidth. 
Its pulse area, proportional to t' l x W' l (Eq. (16)) is equal 
to tt. Therefore, for this pulse, the population transfer 
occurs only toward the resonant level and is complete. 
The corresponding laser sources are not numerous but 
are presently developed. 

Let us emphasize that unchirped femtosecond pulses have 
been considered throughout this paper. Methods for ex- 
ecuting robust, selective and complete transfer of pop- 
ulation between a single level and preselected superpo- 
sitions of levels are presently rapidly developing, both 
theoretically [35] and experimentally [36]. Such tranfers 
are obtained through adiabatic passage with intense fem- 
tosecond pulses, shaped in amplitude and phase. 
We finally mention that the possibilities offered by the 
implementation of a STIRAP (stimulated Raman adia- 
batic passage) process using femtosecond pulses, instead 
of the currently used pulses in the microsecond domain, 
remain to be investigated. Ultimately one would want 
to produce absolute ground state molecules from weakly- 
bound molecules formed after photoassociation and spon- 
taneous radiative decay. Keeping this goal in mind, one 
may want to investigate schemes [37] where a coherent 
train of weak pump-dump pairs of shaped femtosecond 
pulses are used. In that scheme each pair of pump-dump 
pulses drives narrow-band Raman transitions between 
vibrational levels avoiding spontaneous emission losses 
from the intermediate state. 



Appendix A: Complete set of vibrational states from 
the MFGH method 



The MFGH method is based on the Fourier Grid 
Hamiltonian method (FGH) with the introduction of an 
adaptive coordinate, related to the local de Broglie wave- 
length, to represent the interatomic distance R describing 
the vibration of the diatomic molecule in the potential 
V(R). The employed spatial grid has a few points N but 
a large extent L [12]. The Hamiltonian H mo i is repre- 
sented on this grid using a sine expansion rather than 
the usual Fourier expansion, in order to avoid the occur- 
rence of ghost levels [13]. For a single channel problem 
the diagonalization of the Hamiltonian matrix provides 
a complete set of N vibrational wave functions <p v (R) 
(0 < v < N — 1) describing bound levels and discretized 
continuum states normalized to unit on the grid. As dis- 
cussed in Rcf. [17], only a small number of scattering 
wave functions, called 'physical scattering levels', have a 
realistic behavior throughout the grid. The other ones, 
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which have a high probability density at short internu- 
clear distance, ensure the completeness of the set for 
< R < L. The eigenfunctions ip v (R) are orthogonal 
within the box: 



<p Vl (R)(p V2 (R)dR = S VltV2 . 



(Al) 



They satisfy the following closure relations, valid for R < 
L, R' < L: 



N-l 



Y^ <Pv{ R ) x <Pv(K) = ${ R - R ') 



(A2) 



Expanding, at each time t, the wavepackets ^S g (R, t) and 
ty e (R, t) on the stationary vibrational levels f g>v » and 
ip e ,v' of the g and e electronic states (Eq. (13)), we obtain 
the instantaneous amplitude of population, a v « (t) (resp. 
b v /(t)), in the stationary levels gv" (resp. ev'), with 
energies E g<v « (resp. E e y): 



a v " (t) exp 
b v > (t) exp 



E n 



R i 

-I — t — t 



<p gyV »{R)V g (R,t)dR, 



tp e ,v>(R)*e{R,t)dR. 



v=0 



(B4) 



In the present paper, the same spatial grid is used for 
both a 3 S + and 6 3 II electronic states. It contains N — 
511 points with of length L = 1258 ao- The lower and 
excited electronic states possess 48 and 219 bound vibra- 
tional levels. The 'physical scattering levels' describe a 
very small energy domain, less than 0.01 cm -1 , above 
the dissociation limit located at E — 0. The remain- 
ing ones, the 'unphysical scattering levels', cover a large 
energy range up to 35000 cm" . 



Appendix B: Time-dependent study of 
photoexcitation in a diatomic molecule 

Laser-coupled electronic states: 'Wavepacket' 
description 



The evolution of the wavepackets is studied in the ro- 
tating wave approximation (RWA) [38], by introducing 
a frame rotating at the angular frequency lol/2-k, which 
allows one to eliminate rapidly oscillating terms in the 
system of coupled equations (Eq. 8). The new radial 
wave functions corresponding to the lower V g and the 
excited V e potentials, are defined by: 

$ g (R,t) = Xg (R,t)cxp(-lU L t/2), 

* e (R,t) = Xe(R,t)exp(+iui L t/2). (Bl) 

Neglecting the high frequency components ~ ±2wl/27t, 
one obtains the following coupled radial equations: 



d_(* g (R,t) 
m dt\* e (R,t) 



(B2) 



T + V g (R) 

\W eg T + V e (R) 



\w ge 



%(R,t) 

*e(R,t) 



where T is the kinetic energy operator and the potentials 
dressed by the laser frequency V g (R) and V e (R) , are given 
by: 



V g (R) = V g (R) + hw L /2, 
V e (R) = V e (R)-huj L /2. 



2. Laser-coupled vibrational levels: 'Level by Level' 
description 

In the LbyL description, the time-dependent wave func- 
tion \P.(£) is decomposed on the sets g (resp. e m ), with 
g n (resp. e m ) wave functions <p g y (resp. f e . v ') describ- 
ing stationary vibrational levels of the ground and excited 
electronic states. In the interaction picture [38], the ex- 
pression of the wavepackets created by the laser pulse on 
the ground and excited states are given by Eq. (10). 
The time-dependent Schrodinger equation governing the 
time evolution of the ground a v « (£) and excited b v , (t) 
probability amplitudes is equivalent to the system of 
(g-n + e m ) coupled equations: 



a v - 


= 


-\ S ^ ,exp 


ft 


- Eg,v" )t 




X 


0Pg,v" Wge 


Ve,v'J > 




iv 


= 


-j. Yl ^- ex p 

v" eg_ n 


-\i E 9^ 


— E e y)t 




X 


(<Pe,v' \W eg 


■Pg,v» ) • 





(B5) 



In this system, there appear only the off-diagonal matrix 
elements of the coupling W ge (t) = —D ge £of(t)cos[u)Lt]- 
In the RWA approximation, when the high frequency 
[(E e>v > - E g y)/h + u>L] ~ 2lo l and [(E g>v « - E e y)]/h- 
lol ~ —2lol can be neglected, the system reduces to: 



2ft E *« ex p 

4 Yl ^" exp 



ft 



( ( Pg,v"\W ge \tp e y) , 



" es„ 



ft 



(<Pe,v>\Weg\ip g , v 



(B6) 



W ge (R,t) = W eg (R,t) is defined in Eq. (9) and S v iy is 
equal to: 



(B3) 



S v ',v" — ft^L — [E e y — E g y\. 



(B7) 
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If the R- variation of the electric dipole moment D ge (R) 
is neglected, introducing D = D ge (R — > oo), one has: 



(ev'\W eg (R,t)\gv") = ( 5 «"|W 3e (t)|e«') 
= -D£b/(t)<vV) - Vvyfit). 



(B8) 



We solved the differential equations (B6) using the func- 
tion NDSolve of the Mathematica software system. The 
energies and wave functions for the levels v" and v' , as 
well as the overlap integrals (v'\v"), were obtained from 
the MFGH method (see Appendix A). 



3. 'Wavepacket' and 'Level by Level' descriptions 

Using the expansion Eq. (13) of the wavepackets X s (-R, t) 
and Xe(R, t) in terms of the stationary wave functions, 
and accounting for the closure relations satisfied by the 
wave functions <p g y and (p e . v r (Eq. (Al)), one obtains 
a system of 2N first-order differential equations for the 
probability amplitudes of vibrational states a v » (t) and 
b v ' (t) involved in the WP description. This system is very 
similar to the system Eq. (B5) satisfied by the probability 
amplitudes a v „ (t) and b v , (t) in the LbyL description. 
The difference between the two systems arises only from 
the number of involved amplitudes : 2N for the WP 
description and (g n + e rn ) in the LbyL approach. We 
emphasize that the WP description takes automatically 
advantage of the completeness character of the set of 
eigenfunctions provided by the spatial representation of 
the Hamiltonian on a grid. The description of the dy- 
namics does not depend on the choice of the grid pa- 
rameters, provided that a sufficiently wide energy range 
is spanned by the eigenvalues obtained in the diagonal- 
ization. Thus the WP method provides a general non- 
perturbative treatment of the molecule- laser interaction, 
limited to the considered electronic states. It is straight- 
forward to extend the two-states model employed here to 
models with several electronic states. Such multi-surface 
models may become necessary, for example, in studies 
of photoexcitation of vibrational levels belonging to elec- 
tronic states coupled by molecular interactions. 



quency ujl/2t:. The laser is resonant with the transition 
gvo" — > e«Q. The time-dependent wave function vE'(t) 
is explicitly expanded over the diabolic basis made of g n 
wave functions \£ g v") of vibrational levels in the ground 
electronic state and e m wave functions \£ e v 1 ) of levels of 
the excited electronic state: 



*(*) = cx P (-^E g , V0 ;t) ]T A v „(t)\£ g v") 



(CI) 



v"£g 



+ exp(--Eg iVo «t) ^ R v ,(t) exp(~ioj L t)\£ e v). 



v' Ge 



In the RWA approximation, i.e. neglecting the rapidly 
oscillating terms exp(— i2cj^i) (resp. exp(+i2uii,t)), the 
amplitudes A v „ and B_ v , satisfy the system of coupled 
first-order differential equations: 



i A... 



K-,vo'4* + i E <V,«" /(*)£,'. 



i>'£e 



iB v , = K>y B v , + \ Y, ««',«•/(*) 4," 



(C2) 



"" e 9„ 



where Cl v ' tV » are given by Eq. (B8) and where the energy 
differences Ag°^ are defined in Eq. (11). 

The time-dependent Hamiltonian Rdiab(t) is represented 
in the diabatic basis by the following matrix: 





ft-dici 


b(t) 
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9n 






A 9 ,, „ 

v ,v 













A 9 ,, , 

V ,v 















A 9 ,, „ 



o ' *v f ,v" J\P) 



2 ^ ^v',v" J \£) 



A e , , 

v',v 





C3) 



Appendix C: RWA, diabatic and adiabatic basis, 
adiabaticity 

1. RWA at the laser frequency, diabatic basis 

The interaction picture has been used, in the LbyL frame- 
work (Sec. B2), to analyze the dynamics of the vibra- 
tional population transfer. In this approach, the Hamil- 
tonian is non-diagonal, with matrix elements including 
terms exp[±i6 v > tV »t] (Eq. (B6)), with oscillating contribu- 
tions depending on the detuning of the laser with respect 
to the frequency of the gw" — > ev' transition. 
Instead of working in the interaction picture, one may 
transform into a reference frame rotating at the laser fre- 



2. Instantaneous adiabatic basis 

At each time t, the diabatic time-dependent Hamiltonian 
T~Ldiab{t) can be diagonalized, determining the (g n + e m ) 
field-dressed or adiabatic levels \j(t)), with eigenvalues 
ej(t) [24]: 



H diab (t)\j(t))=e j (t)\j(t)). 



(C4) 



These states can be considered as a family of solutions 
of the time- independent Schrodinger equation, with the 
time t as a parameter. The normalization condition is 

(j(t)\j(t)) = 1 and the integral (j(t)\ j(t)\, is thus purely 

imaginary. The phase of each eigenvector can be chosen 
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arbitrarily at each time t, and it is possible to choose the 
phase in such a way that (j(t)\ j(t)\ = [39]. 



The diabatic Hamiltonian can equivalcntly be written in 
the adiabatic basis H a diab(t) = Hdiab(t), with: 

9n+e m 
Ha d iab(t)= J2 ^(t)\.](t))(j(t)\, (C5) 

and the wave function *&(t) can be decomposed on the 
adiabatic basis: 






(C6) 



but this condition is clearly questionable [20] and other 
criteria are given, such as 



h (t)\jm(t) 



j{t) - e m (t) 



< 1, or 



(jk(t)|ffi|j m (t)) 



[ej(t)-e m (t)]2 



<1. 



(Cll) 

In fact, there is no fully general rule allowing to predict 
the validity of the adiabatic approximation [24]. The 
more or less adiabatic character of the evolution can be 
discussed a posteriori, once the evolution of the system 
has been computed, by studying the evolution of the pop- 
ulation in different adiabatic levels. 



The amplitudes of population ej (t) of the instantaneous 
adiabatic levels \j(t)) obey the following system of (g n + 
e m ) coupled equations: 



3n+e„ 



h^±=e j (t)e j (t)-ih J2 «**(*)«*(*)• (C7) 



fe=i 



The coefficient a hk {t) = (j{t)\ k(t)\ = - (j(t) \k(t) 

describes the variation of the adiabatic level \k(t)) in the 
adiabatic basis [24]. With the particular phase conven- 
tion written above [39], the sum over k in Eq. (C7) does 
not include k = j. 
An expression of cxj t k{t) for k =/= j is 

[e fc (t) - eM a jtk (t) = (j(t)\ ^^ !*(*)> ~ ^#W 

(C8) 



01 



01 



4. Case of a two-level system 

In the diabatic basis describing the two-level system 
[\& s , ^ e ], the effective time-dependent Hamiltonian in 
the RWA approximation is (Eq. (C3)): 



Hdiab(t) — 



+0 il(t) 

O(t) 6 



(C12) 



where S denotes the detuning of the laser excitation, and 
f2(i) = 2^e.gf(t) is the instantaneous coupling. The 
probability amplitudes in the diabatic basis satisfy the 
first-order differential system : 



ia(t) = n(t)b(t), 

ib(t) = 5b(t) + il(t)a(t). 



(C13) 



The population transferred to the excited level at time t 

reads 



3. Adiabatic approximation 

In the adiabatic approximation, the second term on the 
r.h.s. of Eq. (C7) is neglected, and the adiabatic ampli- 
tudes evolve as: 



ej(t) = ej(t = 0) exp 



ej(t')dt' 



(C9) 



In this approximation, when the system is at the initial 
time in an instantaneous eigenstate of the Hamiltonian 
at t = 0, let us say \j (t = 0)), i.e. when in Eq. (C6) 
ej (t = 0) = 5jj , the system remains in the instantaneous 
eigenstate that evolves from the initial one, and there is 
no jump toward different instantaneous adiabatic states. 
The validity of the adiabatic approximation has been dis- 
cussed in several papers [24, 40-43] . From Messiah [24] , 
a condition of validity is given by: 



max Efc\ifc (*)b'o(<) 

mm[\(e k (t)-e J0 (t))/h\ 



«1, 



(C10) 



Pdiab(t) = \b(t)\< 



(C14) 



The adiabatic character of the process can be analyzed 
by introducing the instantaneous adiabatic basis. This 
basis, ["J/-, ^+], can be obtained from the unitary matrix 
U(t) that diagonalizes Hdi a b(t) defined in Eq. (CI2) [38]: 



U(t) 
with eigenvalues 



cosfe-^/ 2 -sinfe-^/ 2 



sin f e^' 2 



cosfe^ 2 



1 



f 



£ ± (i) = -<5±-V^ 2 + 4|£W 



9(t) is defined by the relation 

2\Q(t)\ 



tan6»(t) 
and the phase </>(£) by 



6 



o < e < 7T, 



n(t) = |fi(t)|e**<*>, 0<^<2tt. 



(C15) 



(C16) 



(C17) 



(C18) 
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For the excitation by an unchirped Gaussian pulse, one 
can take <j>(t) = 0. The instantaneous adiabatic states 
&-(t) and 1 &+(t) are related to the diabatic ones by: 



!;!§)-»-«>( 



(C19) 



They arc solution of 



n d 



+ 



E_ 

E+ 



m /o -i 

' 2 V 1 



ft</> / cos 6> 
2 I -sin( 



— sint 

— cos( 



(C20) 






In the adiabatic basis, the solution of the time-dcpcndcnt 
Schrodingcr equation can be written as 

*(t) = e"**' [«(*) [*-(*)) + /3(t) |tf+(t)>] . (C21) 

The coupled system for the amplitudes of the instanta- 
neous adiabatic levels is: 

1 



ip = -Vs 2 + 4n 2 p + ^-e (t) a . 



(C22) 



The population in the instantaneous adiabatic levels 
|\&_(t)) and \^> + (t)) may be found as 



p£L(t) = \m\ 2 , 



p 



diab 



(t) = H*)| 2 . 



(C23) 



The second term on the right hand side of Eq. (C22) rep- 
resents the non-adiabatic coupling between the adiabatic 
levels. When in Eq. (C15) <p(t) = 0, the non-adiabatic 
coupling is proportional to 4p The non-adiabatic cou- 
pling can be neglected if and only if 



\(ne) 2 « (e + - E_f 



(C24) 



The evolution is then adiabatic and the instantaneous 
adiabatic levels evolve as: 



*±(t) = *±(t = 0)exp(— - / E±(t')dt') 

frJo 



(C25) 



At t = 0, we assume that a(t = 0) = 1 and j3(t = 0) = 
or, equivalently, that *(t = 0) = *_(i = 0) = * g . In 
the adiabatic approximation, the probability amplitudes 
in the diabatic basis, the so-called 'Rabi oscillations' [38] 
can be simply calculated by using the general prescription 
of Sec. IV C 3a of Rcf. [38]: 

1) at t = 0, project the initial probabilty amplitude, 
defined in the diabatic basis, onto the adiabatic basis, 
using the transformation (Eq. (C19)). 

2) propagate adiabatic states according to Eq. (C25). 



3) project ^(t) on the diabatic basis. 



Appendix D: Resonant excitation from a iV-fold 
degenerate level 

We consider the excitation from a system of N degenerate 
levels \gv"), v" = 1 to N, toward a single level \ev' ) 
(g n = N e m = 1). The relevant Hamiltonian Hdiabit) 
(Eq. (C3)) can be written as: 



T~Ldiab{t) 

(. 



A 





W(t)/2 W(t)/2 





W(t)/2 




W(t)/2 



W{t)/2 

A 



V 



(Dl) 



/ 



where A is a common detuning and we assume that all 
couplings are equal, W(t) — Qiy f(t) for v" = 1 to N. 
The energies of the adiabatic levels are given by: 

ej(t) - 0forl<j<(N-l), (D2) 
1 

2 

1 



e.(t) = e j=N (t) = - A-^A 2 +NW 2 (t) 



v / A 2 +NW 2 (t) 



The instantaneous adiabatic levels \j(t)) are defined by 
their components Vij(t) on the diabatic levels \i), which 
satisfy 



W(t) 



[Vij + 



uw,. 



v N , 



W(t) 



v N 



+i,i 



for 1 < i < N 



\eAt) - A]V N+hj (D3) 



For the degenerate eigenvalues §.j{t), with 1 < j < (N — 
1), orthogonal eigenvectors can be found. For example, 
we may construct the following orthogonal basis 



ViAt) 



V 3 +i,j(t) = 



ViU + 1) 

3 



for l<i<j and 1 < j < (N - 1) 



i + i 

Vi,j{t) = for (j + 2) <i< (N + l). 



(D4) 



V itj = y g , T , for 1 < i < N 



For the two other eigenvalues §^(t), the eigenvectors are 

determined by: 

for j = N and j = N + 1 

and V(jv+i),j = V e . T , with e T V g , T 

and the normaliza tion condit ion: 



2 Ve >T 
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For resonant excitation (A = 0) one has V & 



g-+ 



/2N 



and 



<::.+ 



T 



V2- 



We consider now the resonant excitation from a particu- 
lar sublevel of the TV-fold degenerate level of the ground 
state. If the conditions of adiabatic evolution are satis- 
fied, the evolution of the amplitudes of population A t (t) 
and B_i(t) in the diabatic levels \g,i) and |e, 1) can be 
calculated by using the general method described in Ap- 
pendix C4, because the instantaneous eigenvectors of the 
Hamiltonian in Eq. (Dl) are known. At t = 0, the ini- 
tial wave function \g, 1) is expanded over the adiabatic 
levels \j(t = 0)) (Eq. (C6)) leading to the initial ampli- 
tudes Vi—ij(t = 0). Each adiabatic level evolves then 
according to Eq. (C9). The amplitude of population of 
the diabatic level \i) at time t is equal to 



<*!*(*)) 



]TU 4= i,,(0)U„(i)cxp[/' 
Jo 



-ejit'W] 



Using Eqs. (D5), one obtains the probability amplitudes 
in the ground state A^t) or in the excited state B_i{t): 



Mt) = 

AM = 
Mi) = 



N - 1 

N 

1 



N 



9(t) 



N 



N 
i 



^«*> 



e(t) 



for 2 < i < N, 



where 0(i) = / W(t')dt'. 
Jo 



(D5) 



Notice that time-evolution of a resonantly-excited two- 
level system is fully characterized by 0(t) (Eqs. (17, D5)). 
By contrast, the evolution of our (N + l)-level system is 
governed by the angle y/~NQ(t). The absolute value of 
the probability amplitude of the excited level is reduced 
by the factor y/~N and the population is redistributed 
among N levels of the ground state. For a large number of 
degenerate levels (N — > oo) and in the high field regime 
(0(i) large), the excitation process is blocked and the 
population remains in the initial level. 
Conversely, in the low- field limit 0(i) <C n, the popu- 
lation in the excited level in the (N + l)-level system, 
B_i(t) ~ 5©(i), is weak and identical to the amplitude 
of population in a two-level system resonantly excited by 
the same pulse. Simultaneously, in the ground electronic 
state, there is no change in the amplitude of population 
of the levels A^t) ~ 1 and A r (t) ~ for 2 < i < N. 
Appendix E: Ultrashort pulse train 

The time-dependent electric field describing a coherent 
train of Gaussian pulses is given by [10] 

, jV-l 
£(*) = 2 51 £oCxp(i$ )f{t-qT rep ) 

q=0 

x cxp[i(oj L (t- qT rep ) +qA(f) ce )} 1 (El) 



where N is the number of pulses, $o is a constant phase, 
/ is a Gaussian envelope given by Eq. (3), T rep is the 
pulse repetition time, and A0 ce is the pulse-to-pulse 
carrier-envelope offset phase shift. 



In the energy domain, the spectral distribution of £ (t) is 
obtained from the Fourier transform of Eq. (El) 



<^M = y/(^- w L)cxp 



'../V-l , A , 



WJ rep) 



x exp (i$ 



sin [A/" (A(j) ce - LoT rep ) /2] 
sin [(A0 ce - LoT rep ) /2] 



(E2) 



This distribution consists of a comb of structures 
("teeth") located at equally-spaced frequencies 

Jm — Wljrep i t), 

with spacing f rep = l/T rep ; S is the offset frequency equal 
to S = A(f) ce f rep /(2ir). Each peak m = 0, 1, 2 etc, has a 



maximum of intensity [f(u 



■Ul) 



2 J\f 2 increasing as the 



square of the number of pulses. Here uj m = 2irf m and 
f(uj) denotes the Fourier Transform of f(t). The width 
of each tooth f rep 2Tr/Af decreases with increasing Af. 
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